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' Abstract. We present results from an ISOCAM survey in the two broad band filters LW2 (5-8.5 fim) and LW3 (12-18 fim) of a 

0. 13 square degree coverage of the Serpens Main Cloud Core. A total of 392 sources were detected in the 6.7 pm band and 139 
^ . in the 14.3 [im band to a limiting sensitivity of ~ 2 mjy. We identified 58 Young Stellar Objects (YSOs) with mid-IR excess 

from the single colour index [14.3/6.7], and 8 additional YSOs from the H — KjK — m 61 diagram. Only 32 of these 66 sources 
were previously known to be YSO candidates. Only about 50% of the mid-IR excess sources show excesses in the near-IR 
J - H/H — K diagram. In the 48 square arc minute field covering the central Cloud Core the Class I/Class II number ratio is 
19/18, i.e. about 10 times larger than in other young embedded clusters such as p Ophiuchi or Chamaeleon. The mid-IR fluxes 
• i-H , of the Class I and flat-spectrum sources are found to be on the average larger than those of Class II sources. Stellar luminosities 

■ are estimated for the Class II sample, and its luminosity function is compatible with a coeval population of about 2 Myr which 

follows a three segment power-law IMF. For this age about 20% of the Class lis are found to be young brown dwarf candidates. 
The YSOs are in general strongly clustered, the Class I sources more than the Class II sources, and there is an indication of 
sub-clustering. The sub-clustering of the protostar candidates has a spatial scale of 0.12 pc. These sub-clusters are found along 



the NW-SE oriented ridge and in very good agreement with the location of dense cores traced by millimeter data. The smallest 
clustering scale for the Class II sources is about 0.25 pc, similar to what was found for p Ophiuchi. Our data show evidence 
that star formation in Serpens has proceeded in several phases, and that a "microburst" of star formation has taken place very 
recently, probably within the last 10 5 yrs. 

Key words. Stars: formation - Stars: pre-main-sequence - Stars: luminosity function, mass function - Stars: low-mass, brown 
dwarfs - ISM: Individual Objects: Serpens Cloud Core 
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1. Introduction 

The youngest stellar clusters are found deeply embedded in the 
molecular clouds from which they form. There are several rea- 
sons why very young clusters are particularly interesting for 
statistical studies such as mass functions and spatial distribu- 
tions. Because mass segregation and loss of low mass members 
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due to dynamical evolution h as not had time to develop signif- 
icantly for ages $ 10 yrs JScalol 1998b . the stellar IMF can in 
principle be found for the complete sample, at least for suffi- 
ciently rich clusters. For ages ^ 10 5 yrs the spatial distribution 
should in gross reflect the distribution at birth, which gives im- 
portant input to the studies of cloud fragmentation and cluster 
formation. Only in the youngest regions of low mass star for- 
mation do we find the co-existence of newly born stars and pre- 
stellar clumps, which allows one to compare the mass functions 
of the different evolutionary stages. Low mass stars are more 
luminous when they are young, being either in their protostel- 
lar phase or contracting down the Hayashi track, which permits 
probing lower limiting masses. Severe cloud extinction, how- 
ever, requires sensitive IR mapping at high spatial resolution to 
sample the stellar population of embedded clusters. 

ISOCAM the camera aboard the ISO satellite 
jKessler et alJ Il996l) . provided sens itivity and relatively 
high spatial resolution in the mid-IR (Cesars kv et alJ fT996). 
The two broad band filters LW2 (5-8.5 /mi) and LW3 (12-18 
/mi), designed to avoid the silicate features at 10 and 20 
/mi, were selected to sample the mid-IR Spectral Energy 
Distribution (SED) of Young Stellar Objects (YSOs) in differ- 
ent evolutionary phases. According to the curr ent empirical 
pictur e for the early evolution of low m a ss stars ( Adams et al. 
1987} lLadal 1 19871: lAndre et all 1 19931 lAndre & Montmerld 
19941) . newborn YSOs can be observationally classified 



into 4 main evolutionary classes. Class objects are in the 
deeply embedded main accretion phase (^ 10 4 yrs), and have 
measured circumstellar envelope masses larger than their 
estimated central stellar masses, with overall SEDs resembling 
cold blackbodies and peaking in the far-IR. Class I sources 
(~ 10 5 yrs) are observationally characterised by a broad SED 
with a rising spectral index 1 towards longer wavelengths 
(cir > 0) in the mid-IR. The Class II sources spend some 
10 6 yrs in a phase where most of the circumstellar matter is 
distributed in an optically thick disk, displaying broad SEDs 
with - 1 .6 < ffiR < 0. At ffiR « - 1 .6 the disk turns optically 
thin, and the sources evolve into the (~ 10 7 yrs) Class III stage 
where the mid-IR imprints of a disk eventually disappear. A 
normal stellar photosphere has ffiR = -3. Thus, while Class 
objects are not favourably traced by mid-IR photometry, they 
are expected to be rare. At the other extreme, Class III sources 
cannot generally be distinguished using mid-IR photometry 
since most of them have SEDs similar to normal stellar 
photospheres. But mid-IR photometry from two broad bands, 
as obtained in this study with ISOCAM, is highly efficient 
when it comes to detection and classification of Class I and 
Class II sources. Thus, considering the fact that these latter 
objects constitute the major fraction of the youngest YSOs, 
the ISOCAM surveys provide a better defined sample for 
statistical studies than e.g. n ear-IR survey s for regions with 
very recent star formation (see Prusti 1999). 

This paper presents the results from an ISOCAM survey of 
~ 0.13 square degrees around the Serpens Cloud Core in two 
broad bands centred at 6.7 and 14.3 /mi. This cloud, located at 



b 11 = 5° and F = 3 2° at a distance of 259 ±37 pc JStraizvs et alJ 
ll996t lFestin 1998), comprises a deeply embedded, very young 
cluster with large and spatially inhomogeneous cloud ex- 
tinction, exceeding 50 magnitudes of visual extinction. O nly 
a few sources are detected in the optical JHartiean & Lada 
Il985t iGomez de Castro et alJ fl98"5 iGiovannetti et alJll998l) . 
Serpens contains one of the riche s t known collect i on of 
Class objects JCasali et alJ Il993[ iHurt & B arson vl Il996t 
IWolf-Chase et alJll998t bavis et alJll999l) . an indication that 
this cluster is young and active. On-going star formation is 
also evident from the presence of severa l molecular outflows 
jBallv & Ladal 1 19831 Iwhite et alJ 1 19951 iHuard et alJ Il997t 
Her bst et all 19 97: Davis et al. 1999), p re-stellar cond ensations 
seen as sub-mm sources J Casali et alJ|l993 UMcMullin et alJ 
ll994lTesti & Sargenil998HWilliams & Mversl2000l) . a far-IR 
source (FIRS1) possibl y associated with a non-thermal t riple 
radio continuum source ([Rodriguez et alll989HEiroa & C asali 
1989; Cu rieletalJll993l) . a FUor-like object fapdapp etal 
1996), and jets and knots in the 2.1 /mi H2 line ( Eir oa et al 
19971 lHerbstetalJll997l) . Investigations of the stellar con- 
tent have been made with near-IR surveys llStrom et al 



1975 IChurchwell & Koornnee3ll986t lEiroa & Casalil 1 19921: 



Sogawa et alJl997UGiovannetfi et al. 1998; Kaas 1999a), iden- 
tifying YSOs using different criteria, such as e.g. near-IR ex- 
cesses, association with nebulosities, and variability. 

In this paper we identify new cluster members, character- 
ize the YSOs into Class I, flat-spectrum, and Class II sources, 
estimate a stellar luminosity function for the Class II sample 
and search for a compatible IMF and age, and finally describe 
the spatial distribution of both the protostars and the pre-main 
sequence population in this cluster. 

2. Observations and reductions 

2.1. ISOCAM 

This work presents data from the two ISOCAM star formation 
surveys LNORDH.SURVEY.l and GOLOFSSO.D.SURMC. 
These were surveys within the ISO central programme map- 
ping selected parts of the major nearby star formation regions 
in the two broad band filters LW2 (5-8.5 /mi) and LW3 (12- 
18 //m). Other results on young stellar populations based on 
these surveys comprise the C hamaeleo n I, II and III reg ions 
jNordh et alJll996l:l01ofsson et alJl998llPersi et all200ol) . the 
p Ophiuchi star formatio n region ( Bontemps et al. 2001), the R 
Corona Austra lis regi on (Olofsson et al. 1999), and the L1551 
Taurus region (Galfalk et al. 2004), all observations performed 
basically in the same way. Data reduction methods were gener- 
ally the same for all regions, using a combination of the CIA 2 
(CAM Interactive Analysis) package and own dedicat ed soft - 
ware. For overviews of su rvey results, seelNordh et alJ {l998 ), 
iKaas & Bontemps! i2Q0(i . and lOlofssonl i2Q0(i . 

In the Serpens Core region about 0.13 square degrees 
(deg2) were covered at 6.7 and 14.3 //m in 3 separate, but over- 
lapping main rasters named CE, CW and CS (see Fig. llOl i. In 

2 A joint development by the ESA Astrophysics Division and the 
ISOCAM Consortium led by the ISOCAM PI, C. Cesarsky, Direction 
des Sciences de la Matiere, C.E.A., France 



1 The spectral index is defined as am = dlog(AF A )/(dlogA) and is 
usually calculated between 2.2 /an and 10 or 25 fim. 
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Table 1. Observational parameters, detection statistics and photometric results for each of the 6 ISOCAM rasters CE, CW, CS, 
Dl, D2, D3. See Fig.[lO|and text. 



a(2000) 5(2000) Size 1 T int pfov <n r0 > N del N 6 . 7 N14.3 N bo th lcr 6 . 7 1<j-r 3 
('x') (sec) (") (mJy) (mjy) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


CE 


18 29 48.3 


01 16 04.5 


13x13 


0.28 


3/6 


4x13 


133 


113 


56 


51 


2 


3 


CW 


18 29 06.3 


01 16 01.4 


13x13 


2.1 


6 


7 


152 


139 


43 


38 


0.8 


2 


CS 


18 29 52.5 


01 02 37.9 


13x16.5 


0.28 


6 


4x13 


165 


152 


41 


34 


1.2 


4 


Dl 


18 29 48.7 


01 15 20.5 


1.8x4.6 


2.1 


3 


44 


20 


15 


9 


9 


1 


2 


D2 


18 29 52.3 


01 15 20.8 


1.8x4.6 


2.1 


3 


44 


24 


19 


11 


10 


1 


4.5 


D3 


18 29 57.8 


01 12 55.4 


4.6x1.8 


2.1 


3 


44 


30 


21 


15 


13 


1 


3 


Tot 2 














421 


392 


140 


124 







1 Approximate size since each raster has tagged edges. 

2 The total number is corrected for multiply observed sources in overlapped regions. 



addition, 3 smaller fields within CE (named Dl, D2, and D3) 
were observed at a higher sensitivity. CE covers the well known 
Serpens Cloud Core, CW is a reference region to the west of the 
Core which is substantially opaque in the optica l but w ithout 
appreciable 60 yum IRAS emission (Zhan g et all 1 1988 a), and 
CS is a region directly to the south of the Cloud Core which 
has a peak in the 60 /mi flux. The rasters were always made 
along the right ascension, with about half a frame (90") over- 
lap in a and 24" overlap in 6. For the larger rasters the pixel 
field of view (PFOV) was set to the nominal survey value of 6" 
except for the CE region in LW2, where a PFOV of 3" was se- 
lected because of the risk of otherwise saturating the detector. 
Also in order to avoid saturation, the intrinsic integration time 
was set to 7^ = 0.28 s except for field CW where the nominal 
survey value of T mt — 2.1 s was used. Each position in the sky 
was observed during about 15 s. For the deeper imaging within 
CE the overlap was 72" in a and 5, and each position was ob- 
served during about 92 s. The individual integration time was 
Tint = 2.1 s and a PFOV of 3" was used for better spatial reso- 
lution. See Table[2for an overview. 



2.1.1. Image reduction 



Each raster consists of a cube of frames which is reduced in- 
dividually, and in total 9112 individual frames were analysed. 
The dark current was subtracted using the CAL-G dark from 
the ISOCAM calibration library and, if necessary, further im- 
proved by a sec ond order dark co rrection using a FFT thresh- 
olding method (Sta rck et alJll996ah . Cosmic ray hits were de- 
tected and masked by the multire sol ution median transform 
(MMT) method ( Starc k et al.l ll996b). The transients in the 
time history of each pixel due to the slow response of the 
LW detector were t reated with the IAS inversion method v. 1 .0 
jAbergel et al.lll996albt) . Flat field ima ges were constru cted 
from the observations themselves. See IStarcket al.l (11999) for 
a general description of ISOCAM data processing. 



2.1.2. Point source detection and photometry 

Bright point sources, even well below the saturation level, pro- 
duce strong memory effects which are not entirely taken out by 
the transient correction. Due to the large PFOV which under- 
samples the point spread function (PSF), also remnants of cos- 
mic ray glitches may be mistaken for faint sources. By looking 
at the time history of the candidate source fluxes, however, it is 
easy to distinguish real sources from memory effects, glitches 
or noise. This is done for each individual sky coverage, the re- 
dundancy being 2-6. Both the source detection and the pho- 
tometry was made with the interactive software develop ed by 
our team. See also previous papers on this survey jNordh et al. 
|1996 | lOlofsson et alJll999l IPersi et alJl2000t iBontemps et al. 
l200lh . 

The fluxes of each verified source were measured by aper- 
ture photometry, with nominal aperture radii 1 .5 and 3 pixels, 
for the 6" and 3" PFOV, respectively, and aperture correction 
applied using the empirical PSF. For each redundant observa- 
tion (2-6 overlaps) the flux is the median flux from the N I0 
readouts per sky position. The sky level was estimated from 
the median image to reduce the noise. Uncertainties are esti- 
mated both from the standard deviation around this median (i.e. 
the temporal noise and also a measure of the efficiency of the 
transient correction method) and from the cr of the sky back- 
ground (i.e. the spatial noise). The quoted flux uncertainties 
are these two contributions added in quadrature. The photo- 
metric scatter between the overlaps is estimated and deviating 
measurements are discarded. If a source is affected by the dead 
column (no 24 was disconnected) or is hit by a serious glitch 
or memory effect in one of the redundant overlaps, this one is 
skipped, and the remaining ones are used to estimate the flux. 
If the redundancy is not sufficient (i.e. along raster edges), then 
the fluxes are flagged if affected by the dead column, detector 
edges, glitches, memory effects, close neighbours etc. 

Source positions were calculated from that of the redundant 
images where the source is closest to the centre, in order to di- 
minish the effects of field distortion. The source centre is taken 
to be the peak pixel, and the default ISO pointing was used as 
a first source position value. The ISO position was then com- 
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pare d to near-IR posit i ons of known sourc e s in the literature 
feiroa & Casalil Il992l ISogawa etaDll997t iGiovannetti et alJ 
1998), and ISO positions of bright optical sources to the digital 
sky survey. Bulk offsets found for each raster registration were 
then corrected. The maximum bulk offset found was 5", and 
we estimated an average uncertainty of ±3" in RA and DEC. 
When 2MASS became available we checked the positions of 
61 sources in Table[2] which have 2MASS counterparts and are 
not multiples unresolved by ISOCAM. The median deviation 
between ISOCAM and 2MASS positions is 2.2". One star de- 
viates by as much as 9.2" (ISO-356), and 5 more sources by 
more than 4.5" (these are ISO-29, 207, 272, 357 and 367). 

2.1 .3. Photometric Calibration 

The two broad band filters LW2 and LW3 have defined refer- 
ence wavelengths at 6.7 and 14.3 jum, respectively. The fluxes 
in ADU/s are converted to mJy through the relations 2.32 
and 1.96 ADU/gain/s/mJy for the bands L W2 and LW3, re- 
spectiv ely (from in-orbit latest calibration by Blom maert et alJ 
(2000)). These flux conversions are strictly valid only for 
sources with F v oc y _1 , however, and therefore a small colour 
correction must be applied to the blue sources (cf. Sect. 13.21 
for blue vs. red). This correction is obtained by dividing 
the above conversion factors by 1.05 and 1.02 for LW2 and 
LW3, respectively, yielding the effective factors 2.21 and 1.92 
ADU/gain/s/mJy for blue sources. 

Conversion from flux density to magnitude is de- 
fined as iri6.7 = — 2.5 log(F v (6.7//m)/82.8) and mm^ = 
-2.51og(F v (14.3)/18.9), where F v is given in Jy. The ~ 5 % 
responsivity decrease throughout orbit has not been corrected 
for. 

2.2. Nordic Optical Telescope near-IR imaging 

A 6' x 8' region inside Serp-CE, the area which is usually 
referred to as the Serpens Cloud Core and has been covered 
to a smaller or larger extent by several studies in the near- 
IR feiroa & Casalill992tlSogawa et al.ll997UGiovannetti et alJ 
ll998HKaasll999al) . was mapped deeply in J ( 1 .25 //m), H ( 1 .65 
/im) and K (2.2 fim) in August 1996, only 4 months after the 
ISOCAM observations. See Fig.[H)]for the location of the dif- 
ferent maps. Also, a region to the NW of the JHK field has 
been mapped in the K band, see Fig. [4] but with a total coadded 
integration time from only 30 sec to 1 min. These observations 
were made with the ARcetri Near Infrared CAmera (Arnica) 
at the 2 .56m Nordic Optical Telescope, La Palma. See Kaas 
lll999al) for details on this near-IR dataset. 

2.3. IRAM 30m Telescope Observations 

A 1.3 mm dust continuum mosaic of the Serpens main cloud 
core was taken with the IRAM 30-m telescope equipped 
with the MPIfR 37-channel bolometer array MAMBO-I 
jKrevsa et al.lll998l) during four night observing sessions in 
March 1998. The passband of the MAMBO bolometer array 



has an equivalent width as 70 GHz and is centered at v e ff « 
240 GHz. 

The ~ 19' x 6' mosaic consists of eleven individual on-the- 
fly maps which were obtained in the dual-beam raster mode 
with a scanning velocity of 8"/sec and a spatial sampling of 4" 
in elevation. In this mode, the telescope is scanning continu- 
ously in azimuth along each mapped row while the secondary 
mirror is wobbling in azimuth at frequency of 2 Hz. A wobbler 
throw of 45" or 60" was used. The typical azimuthal size of 
individual maps was 4'. The size of the main beam was mea- 
sured to be ~ 11" (HPBW) on Uranus and other strong point- 
like sources such as quasars. The pointing of the telescope was 
checked every ~ 1 hr using the VLA position o f the strong, 
comp act Class source FIRS1 (good to ~ 0.1" - ICuriel et alJ 
( 1993P); it was found to be accurate to better than ~ 3". The 
zenith atmospheric optical depth, monitored by 'skydips' ev- 
ery ~ 2 hr, was between ~ 0.2 and ~ 0.4. Calibration was 
achieved through on-the-fly mapping a nd on-off observations 
of the primary calibrator Uranus (e.g. iGriffin & Ortor]ll993l 
and references therein). In addition, the Serpens secondary cal- 
ibrator FIRS 1 , which has a 1 .3 mm peak flux density ~ 2.4 Jy in 
an 11 " beam was observed before and after each map. The rela- 
tive calibration was found to be good to within ~ 10% by com- 
paring the individual coverages of each field, while the overall 
absolute calibration uncertainty is estimated to be ~ 20%. 

The dual-beam maps were reduced and combined with 
the IRAM software for bolometer-array data ("NIC"; cf. 
iBroguiere. Neri & Sieyersl \ 1 9951)1 which uses the EKH 
restoration algorithm ( E merson. Klein & Has lam 1979). 



3. ISOCAM results 

3.1. Source statistics, sensitivity and completeness 

Table^gives an overview of the observational parameters and 
the results of the point source photometry for each of the six 
rasters named in col. 1. Columns 2, 3 and 4 give the centre 
position and size of the raster fields. Columns 5, 6 and 7 give 
the unit integration times, the PFOV and the average number 
of readouts per sky position. Column 8 gives the total num- 
ber of source detections. Columns 9 and 10 give the number 
of sources for which photometry was obtained at 6.7 and 14.3 
//m, respectively. Column 1 1 gives the number of sources with 
flux determinations in both of the photometric bands. Columns 
12 and 13 give the measured lcr photometric limits which are 
about 4 and 6 times the read-out-noise for the large and the deep 
fields, respectively. Practically all sources detected at 14.3 fim 
are also detected at 6.7 fim, but there are 15 cases of 14.3 fim 
detections without 6.7 fim counterparts. On the average, only 
30% of the 6.7 fim detections are also detected at 14.3 fim, 
but this value is clearly larger for the CE field. The sensitivity, 
which depends on the PFOV and T; nt as well as the amount of 
nebulosity and bright point sources, is lowered by memory ef- 
fects from bright sources in the CE field and by the presence of 
a nebula in the CS field. For the deep fields (Dl, D2, D3) the 
main limiting factor in terms of sensitivity is the memory effect 
and the fact that the flat fields are based on fewer frames. 
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The deep fields within the CE area are independently re- 
peated observations (observed within the same satellite orbit), 
and therefore permit a check of the photometric repeatability 
for a number of sources. The median values of the individual 
direct scatter between the measured fluxes were found to be: 
about 30% for the 6.7 /mi band measured from 28 sources be- 
tween 4 mJy and 2 Jy, and about 15% for the 14.3 //m band, 
measured from 14 sources between 30 mJy and 5.7 Jy (note dif- 
ferent pixel scale). This corresponds to a median error of 0.13 
dex in the [14.3/6.7] colour index, which is in good agreement 
with the scatter (of sources without mid-IR excesses) around 
the value expected for normal photospheres (see Fig. |5Jl at a 
flux density of 10 mJy in the 14.3 /mi band. 

To estimate how many field stars to expect at a given sen- 
sitivity within the observed fields, the stellar content in a cone 
along the line-of-sight was integrated in steps of 0.02 kpc out 
to a distance of 20 kpc. The Galaxy was represented by an ex- 
ponential disk, a bulge , a halo and a molecular ring, following 
Wainscoat et al . (1992). Absolute magnitudes at 2.2 /mi (M K ) 
and in the 12 /mi IRAS band (M12), local number densities 
and scale heights of the different stellar populations, as well 
as their contribution to the different components of the Galaxy 
were provided by Wain scoat et alJ i 19921) in their model of the 
mid-IR point source sky. Absolute magnitudes at 7 //m, M7, are 
estimated by linear interpolation between Mr and M12. 





log F,(HH (Jy) 



log F„(14.3 M m) (Jy) 





Total sample 




n^T-Th. ftTTi-i i->-' 



log F„(6.7,um) (Jy) 



log F^(14.3^m) (Jy) 



Fig. 1. Histograms of sources at 6.7 /mi (left) and 14.3 /mi 
(right) for the CW field (upper) and the total sample (lower). 
The boldface steps in the lower panels correspond to the popu- 
lation of mid-IR excess sources. The expected number of galac- 
tic sources in each bin, normalised to each field, from mod- 
elling the stellar content in a cone along the line-of-sight is 
inserted for zero cloud extinction (solid lines), Av ~ 5 (dotted 
lines), and Av ~ 10 (dashed lines). 

The upper panels of Fig.n snow the histograms of the 6.7 
//m and 14.3 /mi sources in the CW field, which is free of neb- 
ulosity and practically free of IR excess sources. The bin size 
of 0.2 in log F v corresponds to 0.5 magnitudes. Inserted are 
the model counts of galactic sources at 7 and 12 /mi, scaled 
to the CW field size. The expected source number per bin is 
calculated assuming no cloud extinction (solid line), an aver- 



age extinction of 0.2 magnitudes in the 6.7 /mi band (dotted 
line), and A6.7 = 0.4 magnitudes (dashed line), corresponding 
to roughly Ay ~ 5 and Av ~ 10, respectively (cf. Sects. |4]and|5] 
about extinction). For the CW field an average cloud extinction 
of Av ~ 7 - 8 magnitudes is in agreement with the location of 
the majority of the sources in a DENIS I - J/J - K' dia gram, 
as we ll as an extinction map based on R star counts ( Cambresv 
Il999l) . Comparing the number of observed sources at 14.3 fjm 
with the model expectation at 12 /mi indicates completeness at 
~ 6 mJy or 111143 = 8.7 mag. The observations at 6.7 /mi are 
estimated to be complete to ~ 5 mJy or ni6.7 = 10.6 mag. The 
lower panels of Fig.^show the histograms for the total sample. 
The boldface steps show the contribution of the mid-IR excess 
sources. The dashed line represents the model expectation as- 
suming an average cloud extinction of about 10 magnitudes of 
visual extinction, i.e. A6.7 = 0.41 mag and A14.3 = 0.36 mag 
(cf. Sects. @]and[5li. Thus, the whole sample taken together, and 
allowing for a subtraction of the mid-IR excess sources, indi- 
cates an overall completeness at ~ 6 mJy for the 6.7 /mi band 
and at 8 mJy for the 14.3 //m band (i.e. at ni6.7 = 10.4 mag and 
mi4.3 = 8.5 mag). 




A Normal photoaphe 
Completeness: =8 mJy 



0.1 

K 43 (Jy) 



i 



Fig. 2. The colour index [14.3/6.7] defined as log(F* 4 - 3 /F*- 7 ) 
is shown on the y-axis and the flux F] 4 3 also given in magni- 
tudes mi4 is shown on the x-axis. The SED index between 6.7 
and 14.3 /mi is given on the right hand y-axis. 53 sources with 
substantial IR excess are located above [14.3/6.7] = -0.2 or 
CKj R 14 = —1.6 (dotted line). While 28 of these were previously 
suggested as YSOs, 25 are new YSO candidates (filled circles). 



3.2. Sources with mid-IR excesses 

For the 124 sources with fluxes in both bands we present 
a colour magnitude diagram in Fig. |2] The colour index 
[14.3/6.7], defined as log(F l v 43 /F<?; 7 ), is plotted against F l v 43 . 
This colour index can be converted to the commonly used index 
of the Spectral Energy Distribution (SED): 



7-14 



log(A H F A]4 ) - log(A 7 F Al ) 
\ogA 14 -\ogA 7 



(1) 
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calculated between 6.7 and 14.3 yum and indicated on the right 
hand y-axis. The completeness limit is given with solid lines 
as the combined effect of the completeness in each of the two 
bands. 

The sources tend to separate into two distinct groups, with 
some few intermediate objects. The "red" sources (circles) are 
interpreted as pre-main-sequence (PMS) stars surrounded by 
circumstellar dust. This sample is not believed to be contami- 
nated with galaxies. At the completeness flux level of 8 mJy at 
14.3 fim the expected extracalagtic contamination is about half 
a source within our map coverage and at the level of 3 mJy (the 
faintest source in our sample at 14.3 //m) the contamination is 
below two sources within the mapped region ( Honv 2003). No 
correction for extinction has been made on these numbers, so 
that they should be considered as upper limits to the extragalac- 
tic contamination. 

We have decided to set the division line between mid-IR ex- 
cess objects and "blue" objects (triangles) at [14.3/6.7] = -0.2 
or a 7 lR 14 = -1.6 (dotted line), which corresponds to the classi- 
cal border between Class II and Class III objects (cf Sect|5}- 

Except for a few transition objects, the "blue" group (tri- 
angles) is mainly located at the colour index of normal photo- 
spheres, i.e. [14.3/6.7] = -0.66 or a 1 ]R XA = -3 (dashed line). 
The spread around this value indicates the increasing photo- 
metric uncertainty with decreasing flux, although some of the 
scatter might be real. M giants have intrinsic excesses in the 
colour index [14.3/6.7] of less than 0.1, while late M giants are 
expected to have excesses of up to ~ 0.2 above normal photo- 
spheres. Although late M giants are few, a number of M giants 
are expected in the observed sample of field stars. Thus, the 
slight displacement of this group of sources above the colour of 
normal photospheres could be due to a combination of extinc- 
tion and intrinsic colours of M giants. The effect of extinction 
is small, however. A reddening vector of size corresponding to 
Ak = 3 is indicated in the figure. (More about extinction in 
Sect.il) 

This same dichotomy in colour was found also in 

l 1 1 i 

the Chamaeleo n dark clouds ( iNordh et alJ 119961). in RCrA 
Joiofsson et al.l Il999l) . and in p Qphiuchi jBontemps et alJ 
l200lh . According to the SED index, a 7 lR 14 , the mid-IR ex- 
cess sources are Class II and Class I types of Young Stellar 
Objects (YSOs). In Chamaeleon I about 40% of the mid-IR 
excess sources had been previously classified as Classical T 
Tauri stars (CTTS), and in p Ophiuchi ~ 50% were previ- 
ously known as Class II and Class Is. In Serpens, very few 
sources are optically visible and the IRAS data suffer badly 
from source confusion. The known YSOs are thus a sample 
of sources classified from near-IR excesses, association with 
nebulosity, emission -line stars, ice features, clustering proper- 
ties, and variability (Eiroa & Casali 1992; Hodapp et al. 1996; 
Horrobin et al. l Il997t ISo gawa et al. 1997!; iGiovannetti et alJ 
1998HKaaslll999ah . Previously known YSOs are not found in 
the "blue" group, although Class III sources, young stars with 
marginal or no IR excess, have their locus there. A few objects 
in the transition phase between Class II and Class III are ev- 
ident in the "blue" group. Class III sources without IR excess 
cannot be distinguished from field stars purely on the basis of 



broad band infrared photometry. They can, however, be effi- 
ciently identified through X-ray observations. 

The fact that intrinsic IR excess is so easily distinguished 
from reddening at these wavelengths, enables an unambiguous 
classification of the youngest YSOs on the basis of one single 
colour index [14.3/6.7]. In this way ISOCAM found 53 Class I 
and Class II YSOs, of which only 28 had been previously sug- 
gested as YSO candidates. 

Positions and photometry of the 53 Class I and Class II 
YSOs are listed among other members of the Serpens Core 
cluster in Table |5] The 71 "blue" sources are listed in Table [3] 
Many of these are field stars, but some are likely to be young 
cluster members of Class III type. Fig.^shows that the CW 
region contains as many as 20 (10) sources detected at 6.7 fim 
(14.3 fim) above the expected fiel d star contrib ution according 
to the galactic model by Wai nscoat et alJ 111992). Because there 
are only two IR excess sources in the CW field, of which one is 
very faint, these are practically all belonging to the "blue" sam- 
ple. Since the total region surveyed is about 3 times as large 
as the CW field, perhaps more than 20 "blue" sources from 
Table|3]could be Class Ills belonging to the cluster. 

4. Near-IR and mid-IR excess YSOs 

When it comes to detecting IR excesses, the J band usually 
limits the number of sources in a J - H/H - K diagram because 
of its sensitivity to extinction. For the ISOCAM observations 
the 14.3 fim band is the less sensitive, unless the sources are 
extremely r ed (see limits in Fig.E) . Using the H-K/K-m^j di- 
agram as in Olofss on et alJ II1999). however, we should be able 
to detect more IR excess objects. Fig.[3]shows that practically 
all sources with excesses in the colour index [14.3/6.7] are also 
found to have excesses in their K - m^j index. Approximate 
intrinsic colours of main-sequ ence, giant and su pergiant stars, 
taken from the source table of Wains coat et alJ ll 19921) and in- 
terpolated between 2.2 and 12 fim, are given as boldface curves 
(solid, dotted and dashed, respectively). 

Sources with ISOCAM mid-IR excesses (circles) separate 
well from those without mid-IR excesses (triangles). A red- 
dening vector has been calculated by fitting a line from ori- 
gin through the 4 sources without mid-IR excesses. The slope 
(1.23) is mainly constrained by the star CK 2, which is be- 
lieved to be a background supergiant, see e.g. Casa li & Eiroal 
( 1996). On the basis of a number of background stars in the 
/ - H/H - K diagram, Kaas ( 1999a) found a reddening law of 
the form A,i ~ A. 9 to fit the Serpens data. Extrapolation of this 
law to 6.7 fim would give a slope of 0.83 in the H - K/K - m^j 
diagram, which is in disagreement with all the 4 sour ces with- 
out mid-IR excess. The A,\ ~ A~ 1J law ( Whittet 1988) gives an 
even shallower slope. By assuming, however, the A^ ~ X 9 
law to hold for H and K, and using the empirical slope found 
in Fig. |3] the extinction in the 6.7 fim band is estimated to be 
Ae.7 = 0.41 Ak- This is i n good agre e ment w ith the values 0.37, 
0.35 and 0.47 found bv Ijiang et all J2003h from the ISOGAL 
survey, whose three values depend on the actual form of the 
near-IR extincti on curve applied. 

As shown in Kaas & Bontemps (2000) only about 50% of 
the IR excess sources which display clear excesses in the sin- 
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Table 2. The 77 ISOCAM sources identified as YSOs in the Serpens cluster. 
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317 


29 m 57 s 8 


1°12'52" 


15.30^ 


0.02 


11.47* 


<0.01 


9.45 6 


<0.01 


332 


5 


1039^ 


35 


EC92+EC95 


318 


29'"57;'8 


1°12'37" 






14.22 


0.01 


11.47 


<0.01 


100 / 


9 


105^ 


20 


EC94 


319 


29 m 57 s 7 


1°15'31" 


15.95 


0.03 


12.67 


<0.01 


10.59 


<0.01 


45 


4 






EC93,CK13 


320 


29 m 57'8 


1°12'29" 






15.61 


0.02 


12.91 


<0.01 


23 


7 


19 


5 


EC91 


321 


29""58 s 2 


1°15'22" 


12.97 


<0.01 


10.92 


<0.01 


9.48 


<0.01 


204 


9 


327^ 


5 


CK4,GEL12,EC97 


322 


29'"58 s 3 


1°12'49" 






14.61 


0.01 


12.25 


<0.01 


32 


3 






EC98 


326 


29'"58 ! 7 


1°14'26" 


16.62 


0.06 


13.82 


0.01 


11.56 


<0.01 


202 


39 


335-^ 


8 


EC103 
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ISO 


o^ooo 


<52000 


J 




H 




K 


o~k 


F 6J 


0"6.7 


pU.3 


0"14.3 


Other ID 


# 


18* 




mag 




mag 




mag 




mjy 


mjy 


mjy 


mjy 




327 


29'"58;'9 


1°12'31" 




_ 




_ 


15.41 


0.04 


9.7 


2.0 


8.1 


4.3 


HCE175 


328 


29 m 59 s 2 


1°14'06" 


11.78 


<0.01 


10.36 


<0.01 


9.46 


<0.01 


223 


23 


186 / 


10 


EC105,CK8,GEL13 


330 


29 m 59 I 5 


1°11'59" 










14.37 


0.02 


423 


7 


643 


9 


HB1 


331 


29 m 59 I 7 


1°13'13" 


_ 


_ 


_ 


_ 


15.37 


0.04 


145 


9 


346 


12 




338 


30"W7 


1°13'38" 


13.40 


<0.01 


11.15 


<0.01 


10.07 


<0.01 


16 


4 






EC117,CK6,GEL15 


341 


30 m 01 I l 


1°13'26" 






15.39 


0.02 


12.90 


<0.01 


20 


3 


24 


11 


EC121 


345 


30"'02 s l 


1°14'00" 


17.16 


0.06 


14.74 


0.01 


12.88 


<0.01 


24 


4 


50 


7 


EC125,CK7 


347 


30 m 02 s 8 


1°12'28" 


15.11 


0.02 


11.72 


<0.01 


9.69 


<0.01 


499 c 


12 


658 


13 


EC 129 


348 


30 m 03 s 2 


1°16'17" 


12.14 


<0.01 


10.89 


<0.01 


10.15 


<0.01 


41" 


11 


59'' 


3 


EC135,GGD29 


351 


30 m 04 I 


1°12'38" 


16.51 


0.03 


13.54 


<0.01 


12.13 


<0.01 


8.4 


1.9 


_ 


_ 


EC141 


356 


30 m 05 I 2 


1°12'44" 


16.08 


0.03 


14.51 


0.01 


13.53 


0.01 


2.5 


1.6 






K40,EC152 


357 


30"' 05* 8 


1°06'22" 














10 


1.7 


18 


2.3 




359 


30" 06*4 


1°01'09" 














32 


4 


47 


4 




366 


30 m 07 s 6 


1°12'04" 


12.20 


<0.01 


10.74 


<0.01 


9.92 


<0.01 


53 


7 


64 


5 


STGM8 


367 


3<ro&n 


1°01'42" 














12" 


5 


19" 


4 




370 


30 m 08 J 4 


0°58'48" 














61 


3 


83 


4 


double in I 


379 


30"' 09*3 


1°02'47" 














99 


5 


67 


6 




393 


30'" 1 l; v 2 


1°12'40" 


12.97* 


<0.01 


11.67* 


<0.01 


11.06* 


<0.01 


6.3 


4.1 


11.2 


4.5 




407 


30 m 13 s 9 


1°08'55" 














2.2 


1.3 


4.9 


5.3 





8 H and K band data from the Arnica 1995 map, which is slightly displaced from the 1996 map and therefore includes this object. 
* ISOCAM source is resolved into two sources in the near-IR, and the near-IR fluxes are added. 

Flux measurement affected by proximity to the dead column. If the dead column cannot be avoided by any of the redundant observations, this 
flagging. If source is located on the dead column, a flux measurement is not attempted at all. 

Source close to the detector edge in all redundant observations. 
e Extended source. 

f ISOCAM source is not quite resolved from a bright neighbour. 

8 Galaxy con tamination? For fluxes ~ 3 mjy at 14.3 fim the "red" sample is expected to contain less than two galaxies in our field according to 
lHonvH2003h . 

'" Flux might be affected by memory effects from other sources. 
" Nebulous sky background. 

Empty space means no measurement available in this work, while a hyphen means no detection. 

Identifier acronyms are related to the following references: SVS (Strom et al. 1976); GGD (Gyulbudaghian et al. 1978); HL (Hartigan & Lada 
1985); CK (Churchwell & Koornneef 1986); GEL (Gomez de Castro et al. 1988); CDF (Chavarria-K et al. 1988); EC (Eiroa & Casali 1992); 
SMM (Casali et al. 1993); MMW (McMullin et al. 1994); HHR (Hodapp et al. 1996); HB (Hurt & Barsony 1996); STGM (Sogawa et al. 1997); 
HCE (Horrobin et al. 1997); GCNM (Giovannetti et al. 1998); K (Kaas 1999) 



gle ISOCAM colour index [14.3/6.7] show up as IR excess 
objects in the J - H/H - K diagram. This was demonstrated by 
plotting the same sources in the two diagrams, using a statisti- 
cally significant sample of several hundred sources in Serpens, 
p Ophiuchi, Chamaeleon I, and RCrA. This result implies that 
a sole use of the J - H/H - K diagram will severely underes- 
timate the IR-excess population of YSOs in active star forma- 

I — 1 j 1 

tion regions. Bontemps et al. ( 2001 ) showed that in p Ophiuchi 

only 50% of the 123 Class II sources have near-IR excesses 
large enough to be recognizable in the J - H/H - K diagram. 
Similarly we confirm here specifically for Serpens that the mid- 
IR excess sources cluster along the reddening line when plotted 
in a J - H/H - K diagram, and that only 50% would be recog- 
nized as having IR excesses from JHK data alone. 

The H - K/K - mq diagram is more efficient than the 
J - H/H - K diagram in distinguishing between intrinsic cir- 
cumstellar excess (circles) and reddening (triangles). Thus, by 
sampling the SED a bit more into the mid-IR, the use of two 
colour indices to select sources with intrinsic IR excesses be- 



comes substantially less prone to the influence from cloud ex- 
tinction. Fig. [5] also shows that the mid-IR excesses apparent 
from the [14.3/6.7] index, always appear already at 6.7 fim, 
but only in half of the cases at 2.2 fim. 

In addition to the 53 IR excess sources obtained from Fig|2l 
another eight sources were found to have IR excess from the 
H - K/K - m-i diagram, i.e. displaced by more than 1 cr to the 
right of the reddening band (dotted line). This means we find 
25% more IR excess sources within the centre region than by 
using the [14.3/6.7] index alone. 

Class II sources are most likely Classical T Tauri stars 
(CTTS), but we can not exclude that some of them are weak- 
lined T Tauri stars (WTTS). Among the Class II sources de- 
tected with ISOCAM in Cha maeleon I about 1/3 had been pre- 
viously classified as WTTS JNordh et alJll996l) . For WTTS 
one could attribute the presence of mid-IR excess but a lack 
of near-IR excess to an inner hole in the circumstellar disk 
jMoneti et al.ll 19991) . CTTS, on the other hand, are believed to 
have inner disks, since they show strong Ha emission, which 
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Table 3. The 71 "blue" ISOCAM sources which have [14.3/6.7] < -0.2. Many of these are likely to be field stars, but 
could be young cluster members of Class III type. Colour correction applied. 



ISO q-2000 <5 2 ooo J crj H <r H K cr K F" 1 cr 61 F; A3 o-i 4 3 Other ID 

# 18* mag mag mag mjy mjy mly mly 



2 


28"'41.2 


1°15'15" 








30 


5 


12 


2 




4 


28 m 41.9 


1°14'06" 








35 


5 


6.4 


2.9 




7 


28 m 42.1 


1°14'31" 








26 


3 


7.9 


3.0 




9 


28 m 43.6 


1°16'10" 








19 


3 


8.0 


2.2 




11 


28 m 45.1 


1°10'19" 








10 


2 


3.3 


3.1 




12 


28 m 45.2 


1°16'50" 








57 


8 


11 


2 




16 


28 m 47.0 


1°10'36" 








53 


6 


13 


2 




20 


28 m 48.4 


1°13'14" 








37 


6 


9 


2 




24 


28 m 49.7 


1 = 1 9' 09" 








60 


7 


15 


2 




31 


28 m 52.6 


mo'io" 








5 d 


1 


1.2 


2.8 




42 


28 m 58.3 


l o 10'55" 








20 


2 


3.5 


1.7 


[SCB96] 41 


46 


28 m 59.8 


1°09'58" 








27 


4 


4.8 


2.5 




49 


29 m 0.3 


1°18'17" 








114 


14 


32 / 


4 




61 


29'" 3.4 


1°18'34" 








28 


4 


8 


2 




63 


29'" 3.8 


l o 10'39" 








26 


4 


4.7 


2.5 




65 


29™ 4.6 


1°20'02" 








33 


3 


8 


2 




66 


29"' 4.7 


1°16'15" 








56 


6 


15 


2 




69 


29"' 5.0 


1°22'04" 








28 


3 


7 


2 




74 


29"' 8.1 


1°12'35" 








13 


2 


2.4 


1.2 




80 


29'" 9.6 


1°09'29" 








12 


18 


4.1 d 


2.0 




82 


29™ 10.6 


1°18'30" 








27 


3 


6.3 


1.8 




84 


29 m 10.6 


1°18'06" 








12 


2 


3.8 


2.5 




91 


29 m 12.4 


1°11'20" 








23 


3 


5.3 


1.9 




94 


29 m 13.4 


1°11'59" 








118 


12 


33 


4 




96 


29 m 14.2 


1°21'30" 








31 


4 


8.3 


2.3 




98 


29 m 14.9 


1°21'55" 








79 


13 


19 


4 




102 


29 m 15.8 


1° 9'32" 








128 


22 


34 


5 




104 


29 m 16.9 


1°18'37" 








28 


2 


7.0 


2.0 




112 


29 m 19.9 


1°20'59" 








65 


8 


38 


5 




126 


29™ 25.6 


1°03'50" 








59 


2 


12 


4 




132 


29 m 26.8 


1°11'56" 








21 


2 


5.0 


2.7 




134 


29 m 27.3 


1°12'58" 








15 


2 


1.4 


2.5 


BD+01 3686 


140 


29 m 28.6 


1°10'24" 








51 


2 


15 


2 




144 


29™29.4 


1°14'03" 








49 


7 


12 


2 




145 


29 m 29.7 


1°13'09" 








37 


6 


9.0 


2.5 




153 


29 m 30.9 


1°15'18" 








10 


2 


0.9 


2.1 




166 


29 ra 32.6 


1°09'46" 








20 


2 


7.0 


2.6 




184 


29 m 37.2 


1°09'15" 








125 


3 


37 


6 




189 


29 m 37.7 


rii'i8" 








67 


2 


31 


2 




190 


29 m 37.6 


1°11'30" 








68 


10 


32 


3 




191 


29 m 37.9 


l o 00'36" 








44 c 


2 


10 


3 




208 


29 m 41.0 


1°12'40" 








36 


6 


9 


4 


SVS18 


209 


29™41.1 


1°20'56" 




8.69 


<0.01 


44 


5 


12 


2 




210 


29 m 41.5 


now 








12 


2 


6 


4 




228 


29 m 45.2 


1°18'46" 


10.09 


8.67 


<0.01 


91 


9 


19 


3 


SVS9 


238 


29 m 47.6 


1°05'12" 








1414 


19 


285 


7 


StRS 208 


243 


29 m 48.9 


1 11'43" 11.18 <0.01 


9.53 <0.01 


8.84 


<0.01 


35 


4 


8 


3 


SVS3 


244 


29 m 49.2 


0°58'54" 








20 


3 


T 


3 




255 


29 m 49.9 


1°09'24" 








91 


2 


24 


4 




264 


29 m 51.6 


0°56'21" 








14 


2 


5 


1 




267 


29 m 51.6 


l o 10'32" 








48 


3 


13 


6 


SVS5 


271 


29 m 52.7 


1°04'38" 








24 


6 


4 


6 




290 


29 m 54.3 


1°17'47" 








27 


5 


4 


2 




296 


29 m 55.8 


1=04' 14" 








118 


5 


31 


4 




300 


29 m 56.1 


l o 00'24" 








53 


19 


12 


13 




311 


29 m 57.6 


l o 10'46" 








26 


2 


15 


3 


BD+01 3689B 
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ISO 


C2000 


^2000 


J 




H 




K 




p6J 


0"6.7 


pl4.3 


0"14.3 


Other ID 


# 


18* 




mag 




mag 




mag 




mjy 


mjy 


mJy 


mjy 




324 


29'"58.4 


1°20'26" 














15 


3 


4 


3 


SVS 16 


332 


30"' 0.2 


1°09'47" 














62 


2 


36 


4 




334 


30'" 0.0 


1°21'55" 














160 


10 


40 


3 


SVS 15 


337 


30"' 0.6 


1°15'19" 


19.13 


0.279 


12.61 


0.004 


8.92 


0.001 


340 


8 


95 


6 


CK2,EC118 


339 


30"' 0.8 


1°19'41" 














70 


7 


15 


2 


SVS 13 


342 


30"' 1.6 


0°59'30" 














40 


9 


8 


6 




349 


30"' 3.7 


1° 4'52" 














97 


6 


21 


9 




372 


30'" 8.4 


0°55'29" 














80 


2 


19 


4 




375 


30"' 8.7 


0°55'16" 














98 


3 


28 


5 




378 


30"' 9.0 


1°14'44" 


11.81 


0.003 


9.68 


0.001 


8.67 


0.001 


70 


7 


14 


4 


SVS7 


382 


30"' 9.2 


1°17'53" 














7 


2 


3 


4 




389 


30™ 10.4 


1°19'36" 














18 


2 


5 


4 


BD+01 3693 


401 


30 m 12.1 


1°16'41" 


11.39 


0.003 


9.59 


0.001 


8.64 


0.001 


66 


11 


13 


5 


SVS6 


409 


30 m 14.7 


1°05'24" 














15 


3 


5 


3 




410 


30 m 14.7 


1°06'17" 














22 


2 


7 


3 





* Binary or close companion not quite resolved by ISOCAM. 

c Flux measurement affected by proximity to the dead column. If the dead column cannot be avoided by any of the redundant observations, this 
flagging. If source is located on the dead column, a flux measurement is not attempted at all. 
d Source close to the detector edge in all redundant observations. 
e Extended source. 

Identifier acronyms are related to the following references: SVS (Strom et al. 1976); CK (Churchwell & Koornneef 1986); StRS (Stephenson 
1992); EC (Eiroa & Casali 1992), SCB (Straizys et al. 1996). 



■ Source not detected at 14.3 fim 

®IR excess to [14.3/6.7] ± Q 



^No mid-IR excess 




2 4 6 8 10 

K - m 7 

Fig. 3. The H- K/K- m^j diagram of the Serpens Cloud Core 
with individual error bars for each source, and some H upper 
limits (arrows). The IR excess sources from Fig. |3 are encir- 
cled. The intrinsic colours of main-sequence (bold solid), gi- 
ants (bold dotted), and supergiant (bold dashed) stars cluster 
close to (0,0). A reddening vector (solid line) is estimated on 
the basis of the four sources without [14.3/6.7] excess (trian- 
gles). Eight additional IR excess YSOs (numbered) were found 
in this diagram, i.e. displaced by more than lcr to the right of 
the reddening band (dotted vector). 

is commonly interpreted as a signature of the accretion pro- 
cess onto the surface of the object, but could also arise in stel- 
lar winds. Optical information is scarce in the Serpens Cloud 
Core because of the large cloud extinction, and we do not know 
what fraction of the Class lis are strong Ha emitters. It is well 
known from studies of CTTS locations in the / - H/H - K 



diagram (Lada & Adamslfl99l iMeveret aljfl997h . that many 
CTTS lack detectab le near-IR excesse s; about 40% in Taurus- 
Auriga according to Stro m et al] Jl993h . 

For typical CTTS the near-IR wavelength region is strongly 
dominated by the photospheric emission, and rather large 
amounts of dust hot enough to produce a strong excess at 2 //m 
are therefore needed in order to distinguish intrinsic IR excess 
from the effects of scattering and extinction in the J - H/H - K 
diagram. In a region such as the Serpens Cloud Core, where 7 
of the Class lis in question here were originally prop osed to be 
YSOs due to their association with nebulosity feiroa & C asali 
Il992h . it is likely that scattered light in the J and H bands adds 
to the lack of detectable near-IR excesses as it gives a bluer 
H - K colour index. 

From an ISOCAM sample o f CTTS with mid-IR excesses 
■ I j ,i 

in Chamaeleon I, Comeron et al. ( 2000) found evidence for the 

presence of near-IR excesses to be correlated with luminosity, 
suggesting an incapability of objects with very low tempera- 
tures and luminosities (young brown dwarfs and very low mass 
CTTS) to raise the temperature needed at the inner part of the 
circumstellar disk to produce a de tectable excess at 2 .2 /mi, in 
agreement with model predictions JMever et al.ll997|) . The ma- 
jority of the Serp ens sources are substantially more luminous. 
As pointed out bv lHillenbrand et all Jl998l) . the larger the stel- 
lar radius (i.e. the younger the star is), the more difficult it is 
to separate near-IR excesses from the stellar flux. This could 
perhaps be part of the explanation for the youngest sources 
(e.g. two flat-spectrum sources without near-IR excess). But 
we found all the Class I sources to possess near-IR excesses, 
such that the role of larger radii in these cases seems to be well 
compensated for, probably by their larger disk accretion rates. 
A statistical interpretation based on broad band photometry is 
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Fig. 4. Arnica K band mapping of the region to the NW of the 
JHK field. ISOCAM sources are encircled and their identifi- 
cations given. Sources 159, 202, 216 and EC38 are YSOs with 
[14.3/6.7] colour excess. Sources SVS9 and 209 are without 
such an excess. The remaining sources are not detected at 14.3 
jim. See Sect.lAl 



likely over-simplified, however, and there are many properties 
intrinsic to the star-disk system (such as disk inclination an- 
gle, disk accretion rate, stellar mass and radius) which con- 
tribute to the complexity of the ind ividual YSOs. From the re- 
sults presented here for Serpens, in Bontemps et al. (200l| for 
p Ophiuchi, and in lKaas & Bontempsl J2000I) for the ISOCAM 
star formation surveys in general, it is evident that one should 
be careful in estimating disk frequ encies among YS Os based 
on near-IR excesses only, see also lHaisch et alJ d200lh . Since 
our study samples only IR excess objects, we cannot say any- 
thing about the disk fraction among YSOs in Serpens. 



5. Characterization of the YSO population 

5.1. Large fraction of Class I sources 

We have combined the photometry from ISO CAM with the 
available K band photometry from 1996 3 (Kaas 1999a) and cal- 
culated the two SED indices: a^ R u and a^ R , which are plotted 
against each other in Fig.[5]for the 39 mid-IR excess sources 



and the 6 sources without mid-IR excesses. The index a^ R 4 is 
close to the index originally used to define the three classes I, 
II and III from the shape of their S EP between 2.2 and 10 (or 
25) fxm bv lLada & Wilkin gl d 1 984l) : lLadal J 1 9 87h . Figure0indi- 
cates the loci of these classes in addition to a transitional class 
referred to as "flat spectrum" sources. 

According to t he updated IR spectral classi fication scheme 
jAndre & Montmerlel 19941 Idreene et"al]l 1994b we tentatively 
> 0.3 as Class I sources, those with 

-0.3 < 



define YSOs with aj£ 14 

~ ' 4 < 0.3 as flat-spectrum sources, and objects 



a 



IR 



with -1.6 < a\ R 4 < -0.3 as Class II sources. There is a 



marginal hint of a gap in the distribution of sources along the 
[14.3/6.7] = -0.2 axis at o ^T 14 » 0.3, which wa s seen also 
in the p Ophiuchi sample (Bontemp s et al-lEoOll) . though at 
« 0.5. It is apparent that we cannot distinguish between Class 
III sources and field stars. (The transition sources in the "blue" 




3 This means the central 8'x 6' JHK field and the NW field shown 
in Fig.@ 



Fig. 5. The spectral index a^ R versus a^ R 14 , or equivalently, the 
K - mi IK - mi4 colour diagram, for a subset of the ISOCAM 
sample. The locations of Class I, Class II and Class III sources 
and field stars are indicated. Note the relatively large fraction 
of Class I sources in this diagram. 



group (see Fig. |2j with some mid-IR excesses are all located 
outside the field for which we have K-band observations.) 

IR excess sources without detection at 14.3 //m (see 
Sec. 4.2) can be classified from the approximately linear re- 
lation between a\ R 4 and a^ R found in Fig. [5] In this case 
Class Is are those which have a^ R > 1.2, corresponding to 
K - m-j > 4.8, Class lis are those which have a\ R < -0.25, 
corresponding to K - m-i < 3.2, and flat-spectrum sources are 
those in between. Thus, the total number of sources in each 
category is larger than apparent from Fig. [5] From Fig. [3] we 
have 1 more Class I, 2 more flat-spectrum sources, and 5 more 
Class lis. Three sources (260, 277 and 308) have only K and 
6.7 fim detections, but their very red colours strongly suggest 
membership in the Class I category. 

The number of Class I sources (19) is thus about equal to 
the number of Class II sources (18) in the central part of the 
Serpens cluster. This is very unusual since Class II sources 
outnumber Class I sources by typically 10 to 1 in star for- 
mation regions. If we include the 13 "flat spectrum" sources 
in the Class II group, the Class I/Class II number ratio is still 
as high as 19/31. This is exceptional, also compared to the 
results obtained from ISOCAM survey s in other re gions; in 
Chamaeleon I this number ratio is 5 /42 jKaasll999bl) and in p 
Ophiuchi: 16/123 fBontemps et al. 2001). Such a large popula- 
tion of Class I sources indicates a recent burst of star formation 
in this region and would be in line with t he rich co llection of 
Class objects found in this cluster bv lCasali et alJ 11993) and 
lHurt&Barsonvl ( ll996l) . 

5.2. Reddening effect on the aj^ 14 index? 

While a spectral index between two fixed wavelengths, such 
as 2.2 and 14.3 /mi, is practical and provides a preliminary 
classification for a large number of sources, a truly reddening- 
independent classification is obtainable only when complete 
SEDs are available over the ~ 2 to 100 pm range (or to 1 .3 mm 
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if Class sources are involved). In the following we discuss the 
possibility that sources defined as Class Is from the a^ R 14 index 
could be heavily extinguished Class II or flat-spectrum sources. 

The source CK2 (probably a background supergiant) is 
seen through more than 50 magnitudes of visual extinction 
jKaasl ll999a). and indicates an empirical slope (=1.6) of the 
extinction law in Fig.|5] as do five other sources without mid- 
IR excesses. This reddening line translates to a slope of 0.92 in 
the K-m-i I K—m\4 diagram, which gives the relation A14 = 0.88 
A7 in terms of magnitudes 4 , by using the relation A7 = 0.41A^ 
obtained in Sect. |4] It is remarkable that the group of mid-IR 
excess sources are lined up along the same slope as the one 
outlined by the reddened field stars. In this sense, some of the 
sources classified as Class I objects on the basis of Fig.|5]could 
be interpreted as highly extinguished Class II objects, either 
deeply embedded in the cloud or seen edge-on. Flat-spectrum 
sources and Class lis are generally expected to be optically vis- 
ible objects. In Serpens, however, only one (ISO-Ser-159) of 
the 13 flat-spectrum sources and only three of the 18 Class lis 
are optically visible. This shows the extreme degre e of embed- 
dedne ss, and is in agreement with the C ls O map of White et alJ 
(1995), which suggests that the visual cloud extinction is larger 
than 30 magnitudes everywhere in the NW-SE ridge and may 
reach values as high as Ay ~ 200 in some places. The fact 
that the slopes are so similar suggests that Class Is are either 
merely more embedded Class lis, or that the extinction behaves 
roughly in the same way for Class I envelopes as for the cloud 
in general at these wavelengths. This would be in line with 
Pad gett et all dl999l) who observed three Class I SED YSOs 
with HST/NICMOS and found edge-on disks in all three cases. 

As shown in Fig.[2]the a 7 lR 14 index is rather insensitive to 
extinction. As many as 15 YSOs in our sample have aj R 14 > 0, 
but a 7 lR 14 < 0, among them the well known DEOS. For some 
of these we have more spectral information than the three flux 
points at 2.2, 6.7 and 14.3 fim. The two sources HB1 and 
EC 129 are relatively isolated YSOs both in the near-IR and 
mid-IR images and should not suffer so much from the source 
confusion that sets a limit to the usefulness of IRAS data for 
the o ther sources. We u se HIRES IRAS fluxes and upper limits 
from lHurt & Barsonvl (fl996) and extend our mid-IR SEDs to 
longer wavelengths. Figure|6]shows that both HB1 and EC129 
(HB2) have rising spectra beyond a clear dip in the SED at ~ 
14 pm. These sources have clearly Class I SEDs, but the dip at 
~ 1 4 yum will produce a blu e al R 14 index. 

iHurt & Barsonvl i 19961) divided the flux from their HIRES 
map of IRAS 18272+0114 between the three sources DEOS, 
EC53 and S68N, while ISOCAM resolves six mid-IR excess 
sources (cf Fig.lTD. of which only EC37 has a 7 lR 14 > 0. 

4 The s ame was found in C ha I (unpublished) as well as in p 
Ophiuchi iBontemps et aljEoOll) . This implies that, for all these star 
formation regions, the extinction is slightly larger at 6.7 yum than 
at 14.3 /an, see also lOlofsson et all 1 1999). Hence, we see no min- 
imum in the extinction law in the 4-8 nm interval, as expected for 
standard graphite-silicate mixtures iMathis et alll979llDraine & Led 
Il984l) . But our results agree with the e xtinction law found towards 
the Galactic centr e bv Lutz et"aH fl9961 and the ISOGAL results in 
Jiana et al. (2003) who find A t4 = 0.87 A 7 when using the Rieke & 
Lebofsky law updates ikieke et alll989l) for the near-IR extinction. 



v (Hz) 




\ (Aim) 



Fig. 6. SEDs of HB1 and EC129. The flux points at 25, 60 and 
100 fim are taken from Hurt & Barsony (1996). Note that the 
60 and 100 fim points for EC 129 are upper limits only. 



It is possible that a broad silicate absorption feature and/or 
the presence of H?Q and CO? ices in the 14.3 yum band 
JWhittet et alll996Ude Graauw et alJl 996 ) may cause absorp- 
tion effects that correspond in magnitude to the observed dips 
in the SEDs at ~ 14 /mi, noting that the effect is there for 
the most deeply embedded sources. Recently, ISOCAM-CVF 
spectroscopy of 42 Class I and Class II YSOs in Serpens, p 
Ophiuchi, Chamaeleon and RCrA by Alexand er et alJ (J2003) 
reveals a number of absorption features from 5 to 18 fim. They 
find that the majority of Class I sources fall in their category a, 
having deep absorption features of ices and silicates, and con- 
clude that in the cases of large extinction the continuum spec- 
tral index between 2 and 14 /jm provides a truer value of the 
shape of the underlying continuum than the observed mid-IR 
spectrum . For the sub-sam ple of 20 sources that were studied 
bv [Alexander et all d2003l) in Serpens we found: four Class Is, 
six flat-spectrum sources, four Class II sources, one "blue" 
source, and we did not detect five of their sources. All four 
Class Is, three flat-spectrum sources, and one Class II fall in 
their group a, and two flat-spectrum sources and three Class lis 
in their category b (only SVS2 is of type c). This comparison 
shows overall good agreement, and indicates that we have not 
overestimated the number of Class Is in our study. 

Also, the presence of shocked molecular H? line emission 
in the 6.7 /mi band fe.g. lCabrit et alJl998tlLarsson et alJ2002l) 
included in the measured flux of some of the youngest sources 
could contribute to give a bluer a 7 lR 14 index than expected from 
a dust continuum. 

Recent 2-D mod els of Class I source geometries by 
IWhitnev et all d2003h which include flared disk and bipolar 
cavity produce mid-IR SEDs with very broad dips around 10 
yum, and overall much bluer mid-IR colours than produced by 
simple ID or simplified 2-D models. 

Based on the previous discussion we trust that the a\ R 14 
index gives a more reliable measure of the SED than the a 7 lR 14 
index. We have also looked at the mid-IR fluxes of Class Is ver- 
sus Class lis. If all the 15 YSOs with a\ R u > and a[ R 14 < 
were reddened Class II sources, one would expect them to be 
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Fig. 7. The average flux at 6.7 and 14.3 yum versus a\ R u for a 
subset of the ISOCAM sample. The location of Class I, Class II 
and Class III sources and field stars are indicated. Filled circles 
mark the 15 YSOs with a 2 m u > and a 7 lR 14 < 0. 



on the average fainter than (or at most as bright as) the YSOs 
with ffj^ 14 < 0. We note, however, that there is a slight ten- 
dency that sources with large a\ R u indices also are, on the av- 
erage, the brightest ones in the mid-infrared. Figure Q shows 
the average flux at 6.7 and 14.3 /im versus a^ R 14 , with a red- 
dening vector of Ar = 3 magnitudes inserted. The filled cir- 
cles in the figure mark the 15 sources that have a\ R lA > and 
a[ R u < 0. Excluding the flat-spectrum sources, the median 
fluxes are 0.18 and 0.029 Jy for 15 Class Is and 27 Class lis (in- 
cluding all those from Table[6]that have mid-IR fluxes), respec- 
tively. For comparison, the median mid-IR fluxes of 16 Class Is 
and 7 6 Class lis in p O phiuchi are 1.17 and 0.095 Jy, respec- 
tively JBontemps et aljEoOlh . In Chamaeleon I for 5 Class Is 
and 42 Class lis these numbers are 0.28 and 0.077 Jy, respec- 
tively. While the statistics is low for the Chamaeleon I Class I 
sources, there is an indication of luminosity evolution from 
the Class I to the Class II phase, the case being strongest for 
p Ophiuchi. We conclude that on the basis of this statistical 
luminosity argument, there seems to be an intrinsic difference 
between the Class I and the Class II populations. 

Noting that Class Is are o n the average more lu minous in 
the mid-IR than Class lis (cf. Bontem ps et alJ l2001). we cau- 
tion that it is not entirely excluded that the large number frac- 
tion of Class I sources in Serpens could be a direct effect of a 
lower sensitivity because of a larger distance. There are about 
80 faint YSO candidates in t he Serpens C loud Core below the 
sensitivity limit of ISOCAM JKaasll 999a). 

Although the exact Class I/Class II number ratio might be- 
come subject to modification when future observations at high 
spatial resolution and sensitivity towards longer wavelengths 
are available (SIRTF-Spitzer at 24 fim and FIRST-Herschel at 
90-250 yum), we believe that the main conclusion of an unusu- 
ally large fraction of Class I sources in the Serpens Cloud Core 
will be maintained. 



5.3. The protostar sample in Serpens 

Based on the discussion in the previous paragraph, we have 
arrived at a sample of 20 Class I SED sources which are listed 
in Table|4]with the two SED indices u\ R lA and a\ R ■ These are 
all protostar candidates. Although Class sources normally are 
not expected to be detectable at shorter wavelengths than about 
25 jt/m, we detect two Class I sources, ISO-241 and ISO-308, 
within the + 3" positional uncertainties of S68N and SMM4, 

res pectively. 5 

iMotte & Andri feOOll) found that about 60% of the Class I 
sources in Taurus-Auriga are "true protostars", using the cri- 
terion that the envelope mass to stellar mass fraction is ^ 0.1. 
A similar method has been applied to the Serpens data. But 
first we exclude 3 sources in order to account for the possible 
confusion with Class sources (in the case of EC41) or contri- 
bution of Class sources (in the cases: ISO-241, ISO-308) in 
the sample of 19 Class Is in the region mapped by IRAM. The 
fraction of bona-fide protostars in the Serpens Core is then esti- 
mated to be 9 out of 16 Class Is or 56%, with an uncertainty of 
+ 10% because of source confusion problems with the strong 
clustering. 

Extended, elongated and polarized near-IR emission is pre- 
dicted by models of infalling envelopes which have developed 
cavities owing to bipolar outflows, where near-IR radiation 
from the central object may escape an d scatter off dust in th e 
cavity and the outer en velope, see e.g. IWhitnevetai1lll997l) . 
IPark & Kenvonl d2002l) find that 70% of the bona-fide proto- 
stars in Taurus are extended in the near-IR, and that less than 
10% of the sources that are extended in the near-IR show no 
extension in mm continuum. Our near-IR data has a spatial res- 
olution about 10 times better than the IRAM map, and we have 
included in Table 0] the full width at half maximum (fwhm) 
of a gaussian fit to the source profiles in the K band images. 
The median fwhm of 160 isolated and relatively bright sources 
over the field is 1 ."27 ± 0."10. Sources with a fwhm greater 
than l."6, i.e. above 3cr, are here defined as extended in the 
K-band. For the Class I sources this concerns 8 of 19 sources 
or 42%, highly coinciding with continuum sources except for 
ISO-331 and ISO-260. There are no extended sources in the 
flat-spectrum and Class II samples. 

Substantial K-band variability was found in 8 of the above 
Class I sources, 5 of the fla t-spectrum sources and 5 of the 
Class II sources llKaasll999ah . The brightness variations of the 
Class I sources are given in Table 0] as AK. The median value 
of the amplitude is 0.38 + 0. 27 mag, la rger than the variations 
found for Class Is in Taurus dPark & Kenvon 2002 

The flat-spectrum sources are believed to be in a transition 
phase between the Class I and the Class II stage, but a few of 
these might be found to be true protostars. We have listed them 



5 According to l Andre et alj i2000) the following sources satisfy the 
Class criteria: FIRS1, S68N, SMM3, and SMM4, all of which are 
bright sources in the 1 .3 mm IRAM map in FigJjT] In addition, SMM2 
is a candidate Class 0. ISO-258b, which may be a knot in the extended 
emission (cf. Fig lA.U . is 8" away from FIRS1, but could result from 
scattered light (see Sect.|Aj. SMM2 and SMM3 have no ISOCAM 
detections within the positional uncertainties, that is to a 3<r upper 
limit of 6 and 9 mjy for 6.7 and 14.3 yum, respectively (cf. Table0. 
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Table 4. The 20 Class I sour ces. AK and fwhm refers to mea- 
sured variability 1995-1996 (IKaasI 1 999a) and extendedness in 
the K-band, respectively. K-band magnitudes from 1996 are 
found in Table [2 

ISO Other ID a-„ l * fwhm AK 



# " mag 



270 




2.56 


3.54 


2.31 




331 




2.51 


4.12 


4.57 




330 


HB1 


2.35 


4.26 


1.47 


0.32 


250 


DEOS 


2.06 


3.52 


7.52 


0.92 


276 


GCNM53 


2.02 


3.05 


1.89 




249 


EC37 


1.98 


2.66 


1.25 


0.26 


241 




1.74 


3.26 


1.49 




308 


HCE170/171 1 




2.84 


4.48 


0.64 


253 


EC40 


1.23 


2.76 


1.21 




265 


EC53 


0.98 


2.27 


2.27 


0.59 


259 






1.94 2 






260 






1.88 


1.83 




258a 


EC41 


0.96 


1.81 


1.33 




312 


EC88/89 3 


0.85 


1.73 


9.53 


0.11 


254 


EC38 


0.82 


1.78 


1.00 4 




326 


EC 103 


0.62 


1.27 


1.55 


0.22 


327 


HCE175 


0.53 


1.73 


1.42 




306 


EC80 


0.44 


1.70 


1.32 




277 


EC63 




1.21 


1.26 




313 


GCNM94 




1.19 


1.32 


0.43 



1 both are extended, HCE171=K32 is variable 

2 Ks = 14.48 from 2MASS PSC, 16" south-west of extended IRAS 
18273+0059, possible K extendedness not investigated 

3 both are extended, EC88 is variable 

4 from image serp-nw which has better seeing 

Identifier acronyms as in Table|2| 

Table 5. The 13 flat-spectrum sourc es. AK and fwhm refers 
to measured variability 1995-1996 dKaasl fl999a) and extend- 



edness 


in the K-band, respectively. K-band mag 


;nitudes 


1996 are found in Tabled 










ISO 


Other ID 


2-14 


2-7 

a m 


fwhm 


AK 


# 










mag 


234 


EC26 




0.63 


1.31 




322 


EC98 




0.17 


1.31 




345 


EC 125 


0.26 


0.44 


1.30 




317 


EC92/95 


0.19 


-0.03 


1.35 




159 


IRAS 18269+0116 


0.11 


0.23 


1.00 1 




237 


EC28 


0.07 


0.99 


1.23 




307 


SVS2 


0.06 


-0.06 


1.38 


0.19 


347 


EC 129 


0.06 


0.53 


1.34 




314 


SVS20 (double) 


-0.03 


0.30 


2.74 2 


0.14 


318 


EC94 


-0.05 


0.56 


1.35 




341 


EC121 


-0.14 


0.30 


1.27 


0.21 


294 


EC73 


-0.24 


-0.18 


1.31 


0.66 


320 


EC91 


-0.24 


0.45 


1.33 


0.32 



1 From image serp-nw which has better seeing. 

2 Double source, not extended. 



Identifier acronyms as in Table|2| 



in Table[5]with SED indices, /T-band fwhm and variability am- 
plitude. Both SVS2 and SVS20 have centrosymmetric polar- 
ization patterns, indicating that evacuated bipolar cavitie s sur - 
round them, and both are double sources llHuardetalJll997t) . 
Among the flat-spectrum sources ISO- 159, EC 129, SVS20, 
and ISO-237 (i.e. 30%) are detected as mm continuum sources 
in our IRAM map. 

5.4. The pre-main sequence sample in Serpens 

Based on the ISOCAM YSO sample found in Table [T] we have 
used a combination of criteria to arrive at a tentative sam- 
ple of Class II sources in Serpens. These are listed in Table |6] 
with SED indices, whenever available, and additional criteria 
used to argue for a Class II status. Since our K imaging from 
1996 only covers about 10% of the ISOCAM survey, we have 
also used published K band photometry whenever available to 
calculate the SED indices. Also, we have assumed that if the 
source shows strong mid-IR excess and is optically visible, it is 
more likely to be a Class II than a flat-spectrum or Class I type 
of YSO. As a supporting argument for Class II designation we 
have found for a few mid-IR excess sources without JHK data, 
strong Ha in emission (to be published in a future paper). When 
the 2MASS All-Sky Data Release became available while we 
were finalizing our investigation, we used the near-IR photom- 
etry from their point source catalogue for the Class II source 
candidates where such data was lacking, and found spectral in- 
dices in agreement with our above reasoning. 

The pre-main sequence population also includes the 
Class III type of YSOs, objects which show no IR signatures 
of an optically thick disk. Since the basic criterion for select- 
ing YSOs in our study is IR excess, we do not sample these 
sources. With the existing photometric data alone we cannot 
distinguish Class Ills from field stars, except for a few objects 
in the transition zone between the Class II and Class III phases. 
The clear gap between the two groups formed in the ISOCAM 
colour-magnitude diagram in Fig.|2]demonstrates that there are 
few objects undergoing this transition, i.e. that the transition 
must be rapid. This is evident in all the star formation regions 
surveyed by ISOCAM. In Serpens at least five sources from 
Table [3] belong to this transition group, cf. Fig. [2] and we list 
them as candidate Class Ills in Table [7] None of these are lo- 
cated in the central Cloud Core region covered by Arnica deep 
near-IR imaging. 

6. Luminosity distribution 

Only about 10% of the Serpens area surveyed by ISOCAM 
was covered by deep JHK photometry. The total sample is 
ther efore not as homogeneous as our surveys in p Ophiuchi 
dBontemp s et al.l l 200ll). Ch amaeleon JPersi et alJ l2000l) and 
RCrA llOlofsson et all 19 99). Also, since Serpens is at a larger 
distance the overall sensitivity is lower, and source confusion is 
worse, especially in the very clustered active centre. With our 
deep RIJHK photometric coverage over the whole ISOCAM 
survey, which will be published in a future paper, we expect to 
extend the population of IR excess sources, and reveal YSOs 
which are below the ISO sensitivity. With near-IR photometry 
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Table 6. The 43 ISOCAM YSOs found to be Class II sources or candidates. A hyphen for the SED indices means no detection 
at 14 fim. Ak is the derived extinction in the K-ba nd, and L» is the stellar luminosity estimated as explained in Sect. [6] AK is 
the amplitude of the K-band variability 1995-1996 ( Kaas 1999a). None of the sources for which Arnica K images were available 
was found to be extended in K. 
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-1.12 


2.05 


0.16 


Clear IR excess in the H — K/K - m 7 diagram 
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K40,EC152 




-1.03 


0.70 


0.030 


Clear IR excess in the H — K/K - m 7 diagram 


357 




-1.11 


-1.73 


0.25 


0.096 


Ks =11.23 from2MASS 


359 




-0.92 


-1.20 


0.38 


0.33 


Ks = 10.56 from 2MASS, strong Ha emission (to be publ.) 


366 


STGM8 


-1.08 


-1.29 


0.59 


0.59 




367 




-0.87 


-1.16 


0.34 


0.12 


Ks = 1 1.67 from 2MASS, strong Ha em. 


370 




-0.83 


-1.00 


0.02 


0.54 


Ks = 10.13 from 2MASS, dbl in the optical, both Ha em. (to be publ.) 


379 




-1.95 


-2.24 


0.81 


1.2 


Ks = 8.09 from 2MASS 


393 




-1.44 


-2.27 


0.44 


0.066 


double in K 


407 




-1.24 


-2.13 


0.19 


0.021 


Ks = 12.36 from 2MASS 



Identifier acronyms as in Tab|2| 



for the whole sample we can also improve the luminosity esti- 
mate of the Class II sources, by dereddening t he J-band fluxes, 
analogously to what was done for p Ophiuchi ( Bont emps et all 
2001). At the moment we estimate the luminosities from the 
mid-IR 6.7 /urn band only, because this is the most homoge- 
neous measurement we have. The 6.7 yum luminosity function 
is shown in Fig. [8] for the Class I and flat-spectrum sources 
(lower) and Class II sources (upper). 



Protostars are expected to radiate most of their luminosity 
at longer wavelengths. To derive their bolometric luminosities 
we need to integrate the observed SEDs. As shown in Sect.|S| 
however, the four IRAS sources found in the CE region cor- 
respond each to a number of protostars. For reasons of source 
confusion we have not attempted to make calorimetric lumi- 
nosity estimates of the Class Is in Serp ens, as was done fo r the 
Class I sample in p Ophiuchi, see Bontem ps et all J200ll) who 
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Table 7. Five candidate Class III sources from Table [3] These 
are probably in the transition phase between Class II and III 
having some remaining I R excess. K-band magnitudes for two 
sources are obtained from lKaas et alJdl995h . 



ISO 

# 


Other ID 


U IR 


a 2 -' 4 


a IR 


112 




-1.67 






190 




-1.98 


-1.83 


-1.73 


210 




-2.00 






311 


BD+01 3689B 


-1.73 


-2.21 


-2.54 


332 




-1.72 






Identifier acronyms as in 


Tabled 







found that the typical fraction of the luminosity radiated be- 
tween 6.7 and 14.3 fim for a Class I source is ~ 10%. 

As shown in Fig. [S] the protostars span a range in mid-IR 
luminosities about equal to that of the Class lis. It is clear that 
the protostars are on the average more luminous at 6.7 fim than 
the Class lis. This is expected as a large fraction of the lumi- 
nosity of a protostar is the accretion luminosity. The accretion 
can be partly continuous and partly happen in bursts, so that the 
luminosity of a protostar is far from being a simple function of 
its age and mass. 

Class II sources are found to radiate the bulk of their flux 
in the near-IR. The 6.7 /mi band is expected to be dominated 
by disk emission, e.g. for a simple black-body at T = 3700 
K only about 8% of the luminosity is radiated in this band 
(Bo ntemps et alJl200ll) . Thus we can relate the mid-IR flux to 
the total stellar luminosity by assuming that what we see is 
dust emission from a passive reprocessing circumstellar disk. 
For most Class II sources any contribution to the luminosity 
from an active accretion disk is assumed to be negligible. This 
is supported by measurements of the disk accretion rates of 
CTTS in Taurus, which are fou nd to have a m edian value of 
only 10~ 8 M yr~' (e.g. iGullbring et al.lfl998h . If the dust in 
the disk is distributed roughly as in the ideal case of an infi- 
nite, spatially flat disk, 25 % of the cen tral source luminosity 
is absorbed and reradiated by the dust (Adam s & Shul ll986). 
If the disk is flared, the percentage increases, and if it has an 
inner hole, it decreases. Furthermore, the disk inclination angle 
i (to the line of sight) determines what fraction of the disk lu- 
minosity we observe. In addition, we have to correct for cloud 
extinction. Since many of our Class II sources are not detected 
at 14.3yum we select to us e the 6.7/vm flux only. 

Olof sson et alJ dl999t) calibrated an empirical relation be- 
tween the observed 6.7/vm flux and the stellar luminosity by 
selecting sources with known spectral class and very small ex- 
tinction in RCrA and Chamaeleon I. Correcting this relation to 
the larger distance of Serpens we get here that log(L* /L ) = 
logF^j(Jy) + 0.95. Before we can apply this relation, how- 
ever, we have to correct for cloud extinction. For the 
21 Class lis which have Arnica JHK photometry, we assume 
an intrinsic colour of (J - H)q = 0.85, which is a me dian 
value with small dispersion found for CTTS JStrom et alJl993t 
lMeveretalJll997h . and estimate the extinction in the K-band, 
A K = 0.97x[(J-H)-(J-H)o], applying the A -1 9 extinction law 
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Fig. 8. The 6.7 fim luminosity function for the sample of 
Class Is and flat-spectrum sources (lower) and Class II sources 
(upper). 

for the near-IR found for Serpens JKaaJl999ah . Sources with- 
out J-band photometry are either not detected in J - probably 
owing to extinction - or they are located outside the area ob- 
served in the near-IR. In the first case (ISO-Ser-226, 269, 279) 
we interpolate extinction values of their neighbouring Class II 
sources, and in the second case we use the recently available 
2MASS PSC and a dereddening as above. We do not transform 
from 2MASS Ks to Arnica K, since the difference is less than 
or about equal to the estimated errors in the photometry (0.01- 
0.02 mag). 

All extinction values are listed in Table [6] expressed as 
Ak, but only approximate values are given for ISO-Ser-226, 
269 and 279, since it cannot be known at which depth in the 
cloud these objects reside. In Sect. |4] we found the relation 
A(, j = 0.4 IAk, which is used to correct the 6.7 //m flux. The de- 
rived stellar luminosity for each Class II source is also listed 
in Tabled The uncertainty in the luminosity estimate is a func- 
tion of the uncertainties in: 1) the extinction estimate, 2) the 
distance, 3) the L* vs. F(,j relation itself, of which the last two 
will be systematic errors for the whole sampl e. Compariso n 
with previous luminosity estimates by Eiroa & Casali 
but scaled to the distance we use, shows that for the 13 objects 
that overlap in these two samples, there is a large scatter; the 
fraction L/Lec92 varies from 0.5 to 3.7 but the median value is 
1.1. 

The completeness estimate we found in Sect. 13.11 for de- 
tections in the 6.7 fim band is 6 mJy. Because of the variable 
extinction this cannot be directly translated to a completeness 
for Class II sources in terms of stellar luminosities. We have 
used the average measured extinction of Ak = 1.0 for our 
Class II sample and calculated that the completeness limit is 
at L = O.O8L . 

7. Implications for the IMF in Serpens 

lEiroa & Casali i 1992b first estimated a luminosity function 
(LF) for the young cluster in the Serpens Cloud Core, 
including their 51 identified cluster members. The stellar 
luminosities were obtained from a trapezoidal integration 
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over the detected wavelengths, extrapolated to longer wave- 
lengths and extinction-corrected. This LF showed a pro- 
nounced peak around 1 L and a turnover below 0.2 L Q . 
iGiovannetti et alJ (l998) evaluated synthetic K-band luminos- 
ity functions (KLFs) and found a best fit to their observed KLF 
with two bursts of star formation, one 0.1 Myrs ago and the 
other around 3 Myrs ago, and an underlying Miller & Scalo 
IMF. This was partly based on their finding of a turnover in the 
KLF above K = 14 mag. A later study expanded the number of 
Serpens members and found no evide nce for a tur nover of the 
KLF down to a limit of K — 16 mag ( lKaasll999al) . The weak- 
ness of KLFs is that differential extinction is not taken into ac- 
count, and for Serpens this is especially important as values as 
high as Ay ~ 50 mag have been observed for background stars, 
and Ay ~ 10 mag is expected for the densest regions of the 
NW-SE ridge JWhite et all 19951) . 

By selecting sources in the Class II pha se only, we her e con- 
strain the age spread somewhat (Bontem ps et aljEoOlh . This 
phase is thought to last of the order of a few Myr. This means 
that Class lis can have a similar age spread, but we cannot rule 
out that the Class lis once formed in a burst like we see now 
for the protostars (cf. Sect. [S}. We have here made the sim- 
ple assumption of coeval star formation. In Fig.|9]we show our 
observed LF for the sample of 43 Class II sources (shaded his- 
tograms). The bin width of the histogram is dlogL = 0.325, 
based on a factor of two uncertainty in the luminosity estimate, 
and the histogram has been shifted and the bin size slightly 
varied to check that the presented distribution remains stable. 
The observed LF shows a pronounced peak at L ~ 0.09L Q , but 
this corresponds roughly to the completeness limit we have es- 
timated for this sample (dotted vertical line). The error bars are 
the statistical counting errors (VA0. The ma ny lum inous ob- 
jects of around 1L found bv lEiroa & Casalil (Q992) have dis- 
appeared as the sample has been restricted to Class II sources. 

We have assumed coeval population s and three d iffer- 
ent underlying IMFs: the Salpeter IMF dSalpetej i 19551) . the 
Kro upa, Tout and Gil more three-segment power-law IMF 
jRroupa. Tout. Gilmorel 19931) . and the Scalo three-segment 
power-law IMF JScaldll 9981 hereafter S55, KTG93 and S98, 
respectively. With the IMF in the form: dN / d log M* cc M", 
the Salpeter IMF has one single index: a - -1.35, originally 
determined for the mass interval 0.4 to 10 M but in Fig. |9] 
extrapolated to lower masses for reference. Both the KTG93 
and S98 IMFs have three differe nt values of a for three (dif- 
ferently divided) segments. Bontemps et alJ (2001) found for 
a large sample of Class II objects in p Ophiuchi that the high 
mass end of the mass function was well fitted with the in- 
dex a = -1.7. Also, they found that the IMF starts to flatten 
at Mfiat = O.55M and stays "flat" down to 0.055M o with a 
power-law index afi a , = -0.35. These best fits of the two free 
parameters Mfi a , and aji at result in a mass function close to 
both the KTG93 and S98 IMFs for p Ophiuchi. 

On the basis of the pre - main sequen ce evolutionary models 
o flD'Antona & Mazzitellil(ll994[ll998l) . with the 1998 upgrade 
for low mass stars, we have computed synthetic LFs for each 
IMF and for the four different ages: 0.5, 1, 2 and 3 Myr of a co- 
eval cluster. Each window in Fig.[9]shows the result for one age, 
and the computed LFs are overplotted on the observed LF. The 
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Fig. 9. Synthetic LFs were calculated for 4 different ages 
of a coe val p opulatio n on the basis o f pre-main sequence 
tracks of iD'Antona & Mazzitellil Jl998h and 3 different un- 
derlying IMFs. The number of sources per bin is plotted vs. 
logL, and the synthetic LFs are plotted over the observed 
LF fo r the C lass II sources in Serpens (histogram). Both the 
Scalo (1998) three -segment power-law IMF (solid) and the 
Krou pa. Tout. Gilmorel (119931) three-segment power-law IMF 
(dashed) seem plausible for a coeval population of ~ 2 Myr 
all the way down to the estimated completeness limit at L ~ 
0.08L Q (dotted vertical line). The observed LF is also consis- 
tent with the Salpeter IMF (dashed-dotted) for this age at the 
high luminosity end (it has been extrapolated to low masses for 
reference only). 



peak in the LF which wanders towards lower luminosities with 
increasing age arises because of the deut erium burning phase 
dZinnecker. McCaughrean. Wilkindll993l) . Deuterium burning 
acts like a thermostat and hampers somewhat the contraction 
down the Hayashi track, causing a build-up of sources in a 
given luminosity bin for the case of coeval star formation. The 
peak in the observed LF, however, is approximately at the com- 
pleteness limit of the sample, and we cannot put too much con- 
fidence in it. Nevertheless, it is obvious from the figure that the 
observed LF excludes ages of less than about 1 Myr for this 
Class II population for any of three underlying IMFs, which is 
in agreement with current understanding of YSO evolution. It 
is not possible to distinguish between the S98 IMF (solid line) 
and the KTG93 IMF (dashed), but both of them are plausible 
for a coeval population with an age of about 2 Myr, all the way 
down to the estimated completeness limit. Also the Salpeter 
IMF (dashed-dotted line) is in agreement with the observed LF 
for ages of 2-3 Myr at the high luminosity end. Its extrapolation 
to lower masses is given in the plots just for reference. 

Note, however, that if star formation proceeds in a more 
continuous fashion, the structures in the model LFs will smooth 
out. We have checked the model LFs for the case of continu- 
ous star formation over several intervals from 0.5 to 3 Myrs. 
Although the statistics in the observed LF is poor, it seems we 
can exclude the combination of the above tested IMFs and con- 
tinuous SF over periods of less than 2 Myrs. We also tested the 
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Fig. 10. Spatial distribution of the ISOCAM sources in the Serpens Cloud. Sources with mid-IR excesses (filled circles) are seen 
to cluster strongly in the Cloud Core (field CE), while the CW field contains almost exclusively sources without mid-IR excesses 
(open circles). Indicated are also sources with 6.7 fim detections only (plus signs) and 14.3 fim detections only (crosses). The 
coordinates have epoch J2000. Each of the 3 separate ISOCAM rasters, CE, CW, and CS are indicated, and the small box outlines 
the region where deep JHK imaging is available (Kaas 1999a). 



case of a single star formation burst taking place 2 Myr ago but 
having various durations from 0.2 to 1.0 Myr. The peak in the 
observed LF is not secure enough to discriminate between the 
various burst durations. For the tested IMFs and SF histories, 
all we can conclude is that an age of less than 2 Myr for the 
Class II population seems implausible. 

Corrections for binarity have not been applied, but such a 
correction would increase the populations in the lowest lumi- 
nosity bins. We conclude that with the currently obtained LF 
the Class II mass function in Serp ens is similar to the one in 
p Ophiuchi JBontemps et alJEoOlh . even though flat-spectrum 
sources were included in the p Ophiuchi sample. We note that 
the typical stellar mass found for the Serpens Class II sam- 
ple (i.e. the median) i s 0.17 M G , the same as in p Ophiuchi 
(Bontenip s et alJl200ll) . The total mass of the Class II sample 
in Serpens is 16.3M Q . 

We have also used the iBaraffe et alJ Jl998l) evolution- 
ary tracks and calculated synthetic LFs with the same un- 
derlying IMFs as above and for the ages: 1, 2, 3, and 4 
Myr. This model gives basi cally the same result as that of 
iD'Antona & Mazzitellil lll998l) for the degree of confidence we 



can put in our observed LF, taking the large error bars in the 
histogram into account. 

Adopting the above deduced age of 2 Myr for the Serpens 
Class II population, our study is estimated to be complete to 
about O.17M , but to reach down to 0.03 - 0.04M G . Taking 
0.08M Q as the border between low mass stars and brown 
dwarfs, we find 9 sub-stellar size objects in the Serpens Class II 
sample (those with L* < 0.04L Q in Tab. |6j. These are young 
brown dwarfs (BDs) apparently going through a Class II phase 
just like low mass stars. They all have IR excesses with no ap- 
parent deviations from other Class lis in any aspect, which is 
consistent with the idea that they formed in the same way as 
low-mass stars do. Our sample of BDs is not complete, how- 
ever, and we cannot discuss the frequency of occurrence of 
disks among BDs or set strong constraints on BD formation 
models. 

Based on the age estimate above, we find that ~ 21% of 
the Class II sample studied is made up of free-floating brown 
dwarfs. This is similar to the percentage found in p Ophiuchi 
JBontemps et alJEoOlh . The percentage in terms of mass is, 
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however, only about 3% of the total mass for the whole Class II 
sample. 

8. Spatial distributions 

The spatial distribution of the ISOCAM sources is shown in 
Fig. Altogether with the location of the three main rasters CE, 
CW and CS. The filled circles indicate sources with mid-IR 
excesses (cf. Table |2ji, the open circles are the sources without 
mid-IR excesses (cf. Table [3j, the plus signs are sources with 
6.7 fim fluxes only, and crosses sources with 14.3 fim fluxes 
only. The small box outlines the 8' x 6' field for which deep 
JHK imaging is obtained. The overall spatial distribution of the 
mid-IR excess YSOs shows a very clear concentration along a 
SE-NW oriented ridge within the CE field, mainly coinciding 
with what is known as the Serpens Cloud Core. In the CW field 
there is practically no sign of activity, while the distribution of 
IR-excess YSOs in the CS field is quite scattered. The major- 
ity of the mid-IR excess sources seem to be distributed along a 
curved filament pointing towards a centre of curvature around 
18:29:10, +01:05:00 (J2000). Nothing particular was found at 
this location searching SIMBAD, however, and whether star 
formation along this arc may be externally triggered remains 
speculative. The IRAS Point Source Catalogue contains six 
sources in the area surveyed by ISOCAM (cf. AppendixlAlfor 
details). 




Fig. 12. A K ban d image of the Serpens Cloud Core (6' 
x 8') from iKaasI {l999a) with the Class Is (circles), flat- 
spectrum sources (squares) and Class lis (hexagonals) marked. 
The "blue" ISOCAM sources are marked with triangles. The 
ridge along SE-NW is well outlined by the Class I and flat- 
spectrum sources. 



8.1. Clustering scales 

The concentration of young sources found along the NW-SE 
ridge, corresponds well to the pronounced density enhance- 
ment seen in the IRAM 1300 fim map in Fig. \H\ It is evi- 
dent that the Class I and the flat-spectrum sources are spatially 



more confined to the NW-SE ridge than the Class II sources, 
which have a more scattered distribution. In particular, the 
protostars seem to form a set of sub-clusters lined up along 
the ridge. These sub-clusters are in good positional agreement 
wi th the kinematica lly separated N2H + cores A,B,C,D found 
bv lTesti et aD J200(A . especially with B,C, and D. The SVS4 
sub-cluster lies between cores A and B, however, and is also 
displaced to the west of the SMM4 mm continuum peak. 
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Fig. 13. The distribution of source separations for Class I and 
flat-spectrum objects (circles) and Class II sources (dots). The 
number of separations between sources (per bin) is plotted on 
the y-axis as a function of separation in arc minutes on the x- 
axis. The chosen bin size was 6 times the minimum separation, 
translating to 46.1" and 66.1" for the protostar and Class II 
source samples, respectively. 

The spatial distribution of all the classes of YSOs within 
the central 6' x 8' of the Cloud Core can be seen overlaid 
on a K band image in Fig. El To quantify the scale of sub- 
clustering of the two populations: protostar candidates (Class I 
and flat-spectrum sources) and Class II sources, we show for 
each population the distribution of separations between the 
sources. Figure [T3l shows the result as the number of pairs vs. 
separation. The bin size was chosen to be 6 times the mini- 
mum separation. The difference between the two populations is 
clear. The protostar candidates (circles) have a minimum scale 
of clustering of about 0.12 pc (95"), a secondary peak at 0.36 
pc, and practically no scattered distribution, while the Class II 
sources have a minimum clustering scale of about 0.25 pc, a 
secondary peak at 0.59 pc and a much more distributed popula- 
tion, in agreement with the maps. The clustering scale of 0.12 
pc for the protostars agrees well with the range of ra dii of the 
N 2 H + cores: 0.055 to 0.115 pc (Olm i & Testa Eo02). and the 
second peak at about 0.4 pc reflects well the size scale of the 
most active region in the Serpens Core. 

The correspondence between the high density gas and 
the compactness of the protostar clusters indicates that these 
sources are found close to their birth place. Assuming a typi- 
cal velocity dispersion in the range from 0. 1 to 1 km s , YSOs 
younger than 10 5 yrs can be expected to have moved at most 
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0. 1 pc from their birth place. The typical ages of Class I sources 
are estim ated to be of this o rder in other regions such as p 
Ophiuchi jGreene et all 199 4) and Taurus jKenvon et alJl990t> . 
Our results indicate that the protostar candidates in each cluster 
must have formed from the same core at about the same time, 
which may put constraints on models of cloud fragmentation 
and core formation within clumps. Each sub-cluster contains 
between 6 and 12 protostar candidates, and this is just a lower 
limit because of the spatial resolution of ISOCAM. Adopting 
the Serpens distance of 260 pc yields a protostar surface den- 
sity in the sub-clusters in the range from 500-1100 pc~ 2 . To 
our knowledge, such a high spatial density of protostars has 
not been found in any other nearby star formation region, al- 
though we note that the Class surface de nsity in Orion-OM C3 
is comparable to that in the Serpens Core JChini et al J 19971) . 

The sub-clustering in the observed 2D projection of the 
cluster does not necessarily correspond to the true spatial dis- 
tribution of the objects, but it seems highly unlikely that these 
concentrations are due to elongated cloud structures seen end- 
on. Sub-cl ustering was also foun d for the Class II sources in p 
Ophiuchi (Bonte mps et alJl200lh . with a similar size scale as 
for the Class II population in Serpens. It seems likely that the 
strong clustering of Serpens protostars will evolve into looser 
clusters of Class lis over a few 10 6 yrs. The expansion of the 
clustering scales from 0.12 to 0.25 pc and the assumed ages 
of these two populations suggest that the velocity dispersion of 
the young stars is of the order of only 0.05 km s -1 . 

9. Star formation rate and efficiency 

Assuming that a recent burst of star formation took place ~ 2x 
10 5 yrs ago, producing the current population of protostars, we 
can derive a rough estimate of the star formation rate (SFR) and 
the star formation efficiency (SFE) in this burst. We use our 
sample of protostar candidates, the flat-spectrum and Class I 
sources found clustered along the NW-SE oriented ridge. We 
are not able to estimate their masses, but we hypothesize that 
they follow the same IMF as the Class II sample. Comparing 
the total mass of the Class II sample and correcting for the 
number of sources, we estimate a total mass of 12.1 M for 
the protostar population. If these sources were formed grad- 
ually over 2x 10 5 yrs, the SFR in this microburst would be 
6.1 x 10~ 5 M G /yr, or for typical masses of 0.17 M , one new- 
born star every ~ 2800 yr. This SFR is significantly higher than 
the one found for the whole YSO population in p Ophiuchi 
(Bontemps et a l.l200ll) . 

Comparing the protostellar masses with the gaseous mass 
in these cores will give us the star formation effic iency, defi ned 
as SF E = M star /(M star + M gas ). According to lOlmi & T esti 
ll2002l) the two sub-clumps they name NW and SE, which 
contain our four sub-clusters of protostars, are virialised with 
masses of 60 M each. For our protostar candidates this yields 
a SFE of about 9 %. This is the local SFE for the sub-clusters 
along the NW-SE ridge. The global SFE calculated for the 
Class lis and protostars over the entire cluster (i.e. 28.8 M of 
stellar mass) is much lower. In the li terature we find estimates 
of the total gas mass from 300 M Q (Mc Mullinetai]l2000h to 
1500 M (White et al. 1995) for surveyed areas of 5 to 10 arc 



minutes, which gives only upper limits on the SFE of the order 
of 2-9%. 

The above estimates are based on the hypothesis that the 
protostars follow the same IMF as the one found for the Class II 
sources. We have not attempted to derive the protostar mass 
function. Future high resolution far-IR observations (e.g. ESA's 
Herschel Space Observatory) would be needed to yield bolo- 
metric luminosities of these clustered protostars, and with some 
assumptions on the accretion rates one could estimate their 
masses. 

10. Summary and conclusions 

We have used ISOCAM to survey 0.13 square degrees of the 
Serpens Cloud Core in two broad bands centred at 6.7 and 14.3 
jum. In combination with our ground based deep JHK imaging 
of the 48 square arcminute central region as well as additional 
K band imaging of about 30 square arcminutes to the NW, and 
a 1.3 mm IRAM map of the central dense filament, we have 
investigated the mid-IR properties of the young stellar popula- 
tion. The following results are found: 

- The number of point sources with reliable flux measure- 
ments are 392 at 6.7 pm and 139 at 14.3 pm. Of these, 124 
are detected in both bands. 

- On the basis of one single colour index, [14.3/6.7], we 
found that 53 of the 124 objects possess strong mid-IR ex- 
cesses. Only 28 of these were previously suggested YSO 
candidates. The large scale spatial distribution of these 
mid-IR excess sources is strongly concentrated towards the 
Cloud Core, where it is elongated along NW-SE. 

- The near-IR J - H/H - K diagram is found to have an effi- 
ciency of less than 50% in detecting IR excess sources. This 
efficieny is comparable to that one foun d in star formation 
regio ns in general from ISOCAM d ata JKaas & Bontemps 
2000) and p Ophiuchi in particular JBontempseta!Tl200ll) 
and should be kept in mind when interpreting JHK based 
data in terms of disk fractions. 

- The H-K/K-mj diagram separates well intrinsic IR excess 
from the effects of reddening. From this diagram we were 
able to increase the number of mid-IR excess sources to 70. 
This means a fractional increase in the investigated region 
by 25%. 

- Combination of near-IR and mid-IR photometry for red- 
dened stars without IR excesses enables us to estimate 
the extinction at 6.7 and 14.3 pm relative to that in the 
K band in the Serpens direction. We find A7 = 0.41 Ak 
and A14 = 0.36 A^. Our results agree with the extinction 
law m easured towards the Galactic centre by iLutz et all 
dl996 | ). as well as with the results from the ISOGAL sur- 
vey dJiang et alJ l2003h . 

- Classification of the Serpens YSOs in terms of the SED 
indices gives 20 Class I sources, 13 flat-spectrum sources, 
and 43 Class II sources. The number of Class I sources ap- 
pears to be exceptionally large along the NW-SE oriented 
ridge, and the number fraction Class I/Class II is almost 10 
times higher than normal, an indication that this part of the 
cluster is extremely young and active. 
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- The mid-IR luminosities of the Class I sources are on the 
average larger than those of the Class II sources. Since 
Class II sources are expected to have SEDs which peak ap- 
proximately in the near and mid-IR, in contrast to Class Is 
which radiate most of their luminosity in the far-IR, we 
conclude that there is a weak indication of luminosity evo- 
lution in the Serpens Cloud Core. 

- We have estimated extinction and stellar luminosities for 
the 43 Class II sources found in our survey. The Class II lu- 
minosity function is found to be compatible with co-eval 
formation about 2 Myrs ago and an underlyin g IMF of 
the three-segment power-law type dKroupa. Tout. Gilmord 
ll993tLScaloll993l. similar to the mass function found in p 
Ophiuchi (Bontemps et al 12001b . With this assumption on 
age, every fifth Class II is a young brown dwarf. 

- Except for one case, the Class I sources are exclusively 
found in sub-clusters of sizes ~ 0.12 pc distributed along 
the NW-SE oriented ridge. The sub-clusters also contain 
several flat-spectrum sources. In total, each core has formed 
between 6 and 12 protostars (lower limit) within a very 
short time. The spatial distribution of the Class II sources, 
on the other hand, is in general much more dispersed. 

- On the assumption that the protostar candidates follow 
roughly the same IMF as the Class lis, we derive a SFR of 
6.1 x 10~ 5 M /yr and a local SFE of 9% in the recent mi- 
croburst of star formation forming the dense sub-clusters of 
protostars. 

The results presented in this study show evidence that the 
sub-clusters in the central part of the Serpens Cloud Core 
were formed by a recent microburst of star formation. The ex- 
treme youth of this burst, deduced from the compact clusters 
of protostar candidates, is supported by independent investiga- 
tions, such as the rich colle ction of Class so urces found by 
ICasali et all i 19931) and lHurt & BarsonvMl996l) in the same re- 
gions. In addition to the clustered protostar population, we also 
find a more distributed population of Class II sources, for which 
we have deduced an age more or less typical of Class lis found 
in other regions (2 Myr). In addition to these YSO generations, 
there is probably also a population of Class III sources, which 
is undistinguishable from the field star population in our study, 
but which should be loo ked for with prop er search tools, such 
as e.g. X-ray mapping JGrosso et all2000l) . 

ICasali et"a"fl dl993l) suggested that their submm sources 
without near-IR counterparts represented a second phase of ac- 
tive star formation in Serpens. Here we show the co-existence 
of sources in the various evolutionary stages from Class II and 
flat-spectrum to Class I and Class sources. While we find an 
age estimate of 2 Myr for the Class II population, it is highly 
unlikely that the clustered Class I and flat-spectrum sources 
are older than a few 10 s yrs. Our results therefore support the 
conclusions of Cas ali et al.l (11993) that star formation has pro- 
ceeded in several phases in Serpens. 
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Appendix A: Notes on some individual sources 

A.1. IRAS 18269+0116 (ISO-159, 160) 

At the location of IRAS 18269+0116 we find ISO-159 and 
ISO-160, of which the latter is not detected at 14.3 yum. A 
/f-band image shows that ISO-159 has a bow-shaped nebu- 
losity 7" to WNW (at a level of about 3-4 cr, see Fig. |4j. 
The SED index classifies ISO-159 as a flat-spectrum source, 
while ISO-160 is either an extremely extinguished background 
source (analogous to CK2) or most likely a Class II object (cf 
Fig.0. 

A.2. IRAS 18271+0102 (ISO-221) 

Identified as a YSO bv lClarkl Jl99ll) on the basis of its steeply 
rising SED towards longer wavelengths. With ISOCAM we see 
two very bright point sources (221 and 238), of which only 
ISO-221 has IR-excess. 

A.3. EC21 (ISO-231) 

This source was suggested as a YSO on the basis o f its lo - 
cation in the K vs H - K diagram by lEiroa & Casali| (^992). 
Deeper near-IR imaging dete cted the sourc e also in the J band, 
but found no near-IR excess llKaaslll999al) . It is located in the 
very dense and active S68N region, but whether it has intrinsic 
IR-excess at 7 yum or is only reddened (in that case, by 3 magni- 
tudes in the K band, which translates to Ay ~ 30 magnitudes), 
is not clear. 

A.4. EC23 (ISO-232) 

Proposed as a YSO on the basis of its H - K colour and 3.08 
/i/m ice feature llEiroa&Casalilll992l) . but not found to have 
near-IR excess in a deeper JHK study jKaaslfl9"9 9a). nor mid- 
IR excess in this study. 
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A.5. IRAS 18272+0114 (DEOS, EC37, EC38, EC53, 
234, 237, 241) 

IRAS 18272+0114 is resolved into a cluster of protostar can- 
didate s: the Deeply Embedded Outburst Source jHodapp et alJ 
Il996l) . EC37, EC38 and EC53 dEiroa & Casalill992l) . and three 
new ISOCAM sources (234, 237, 241). ISO-241 is within the 
positional uncertai nties associated with S68N, a radio source 
detected by iMcMullin et alJ ll 19941) and found to satisfy the 
Class cr iteria. We use the po sition at its 450 /mi emis- 
sion peak dWolf-Chase et al.ll998l) . ISO-241 is a Clas s I, how- 
ever, f rom the shape of its SED from 2 to 14 /mi. Hodapp 
i 19991) denotes the ZT-band source SMM9-IR, since it is lo- 
cated close t o the SMM9 positions quoted in the literature 
JWhite et alJ Il995t IWolf-Chase et alJll998l iTesti & Sargen3 
1998; Davis et alJI 19991) . and shows that it has a very faint K 
companion about 2" to the NE. 



A.6. EC41 (ISO-258a) and FIRS1 (ISO-258b) 

IRAS 18273+0113, a lso called F I RS1 a nd SMM1 
jRodriguez et al.1 Il989t lEiroa & Casalil Il989t ICasali et alJ 
1993), is associated with 5 mid-IR excess sources (EC41, 253, 
270, 276, 277) plus one 14.3 /mi source (258b). 

ISO-Ser-258a is i n agreement with the near-IR source 
positions of EC41 feiroa & Casalil 1 19921) and GCNM23 
(Gio vannetti et al.lll998l) . and these are most likely the same 
source, a Class I YSO according to the Q'[ R 14 index. The 
ISOCAM source appeared slightly extended to wards the SE 
both a t 6.7 and 14.3 fjm in most of our data. lLarsson et alJ 
(2000) found evidence for two separate sources, however, trac- 
ing pixel by pixel the 5 to 15 /mi SED from CVF observations. 
This prompted us to take a more careful look, and we identi- 
fied two clearly separated sources at 14.3 /mi (rather than one 
extended source) in one of the deeper (and higher resolution) 
survey maps, D2 (cf. Fig. lA.H . The new source (called 258b), 
is located about 13" to the SE of EC41. Since the radio po- 
sition at RA(2000) = 18 29 49.73, DEC(2000) = 01 15 20.8 
dCurielet al.lll993l) is 8.6" to the west of 258b we find it un- 
likely that this source is the mid-IR counterpart of the far-IR 
source FIRS1, but it could be scattered light we see in the mid- 
IR. It is not entirely clear if this source is just a knot of the 
extended emission which appears to begin somewhat south of 
DEOS and is seen in both filters in Fig. IA.1I The northern- 
most part of this extended emission is, however, disturbed by 
a memory feature from the strong DEOS source in the raster 
scanning. 



A.7. IRAS 18273+0059 (ISO-259) 

A prominent extended emission seen in the ISO CAM image s 
corresponds in position to IRAS 18273+0059 Jciarklll99ll) . 
With ISOCAM, however, we see no apparent bright point 
source embedded in the nebulosity, but a moderately faint ob- 
ject is detected at 6.7 /mi only. On optical plates the SE part 
of this extended emission is recognised as an optical reflection 
nebula. 
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Fig. A.l. The deep map D2 at 6.7 /mi (left) and 14.3 /mi (right). 
North up and East left. 



A.8. IRAS 18274+0112 

IRAS 18274+01 12 is also composed of multiple mid-IR excess 
sources centred on SVS20 (cf. FigLL2J. Hig her spa tial 25 /mi 
resolution was obtained bv iHurt & B arson vl ILU?96), and their 
sources PS1 and PS2 (EC 129) are among the bright mid-IR 
excess sources resolved with ISOCAM. 



A.9. EC59 (ISO-272) 

Included as a YSO candidate solely on the b asis of its location 
in the K vs H — K diagram ( Eiroa & Ca sali 1992), and no near- 
IR excess found in a deeper .THK study JKaasll99 9a). Since we 
have no indication of mid-IR excess its status remains incon- 
clusive. 



A.10. EC69 ISO-289 

This source (also called CK10) was suggested as a YSO 

on the basis of its location in the K vs H - K diagram 

■ ■ ■ 

(Eiroa & Casali 1992), and while near-IR excess was found by 



Sogawa et al. ( 1997), no near-IR excess was found in a deeper 
JHK study l lKaaslll999al) . The J-band magnitude has varied 
though: 15.7 in 1989, 14. 96 in May 1992, 15.6 in Nov 1992 
and 16.15 in Aug 1996 feiroa & Casalil 1 19921 ISogawa et all 
ll997tlGiovannetti et alJll998UKaaslll999al) . We conclude that 
it is a YSO, and with an a^ R 7 = -1.42 it belongs in the Class II 
group. 
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A. 11. EC95 & EC92 (ISO-317) 

ISO-317 is located at the position of EC92 a ccording to the 
coordinates given by Eir pa & Cas aljl 1 119921) . but the posi- 
tional uncertainty is roughly +3". Also, the spatial resolution 
of ISOCAM is not sufficient to resolve the two neighbours 
EC92 and EC95, and it is likely that both are included in the 
ISOCAM fluxes. The mid-IR SED suggests a flat-spectrum 
source. Slightly closer to the position of EC95, but also with 
some positional uncertainty, Preibisch ( 1998) found an ex- 
tremely strong X-ray source (Ser-X3). 



A.12. ISO-331 

This bright mid-IR source is detectable in K as a nebulous spot 
only. Probably the K flux is dominated by scattered light, and 
it is therefore not entirely correct to place ISO-331 in colour 
diagrams together with continuum sources. Nevertheless, the 
K measurement gives a lower limit to the colour index of the 
object, which is classified as a Class I sou rce. Its mid-IR po - 
sition is about 25"to the NNW of SMM2 JCasali et alJI 199 A 
and coincides roug hly with the VLA 3.6 cm continuum source 
#11 found bvlBontempsI Jl996l) . 



A. 13. CK2(ISO-337) 

According to ISOCAM data the source CK2, which is ex- 
tremely red in the near-IR, has no mid-IR excess and can be 
interprete d as a background source , consis t ent with the sugges - 
tions of IChurchwell & Koornneed Jl986 | ). IChiar et alJ |l994), 
ICasali & Eh~oaIir996l) and lSogawaetalT i 1997I) . The extinction 
towards CK2 is high, Ay > 50 mag. 
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Fig. 11. Dust continuum mosaic (contours and greyscale) of the Serpens main core taken at 1 .3 mm with the IRAM 30m telescope 
and the MPIfR 37-channel bolometer array (MAMBO). Contour levels: 50. to 200 by 50, 300 to 1000. by 100. mJy/1 l"-beam; 
rms noise level: ~ 17 mJy/ll"-beam. The Class sources are indicated with names and arrows. The location of Class I sources 
(red circles), fiat-spectrum sources (red crosses), and Class II sources (blue filled squares) demonstrates the strong clustering of 
protostar candidates along the dense filament. 



